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A semi-quantitative cluster approach is developed to describe the charge transfer (CT) electron- 
hole excitations in insulating cuprates in a rather wide energy range up to 10 ^ 15 eV. It generalizes 
the Zhang-Ng (ZN) model of CT excitons by considering the complete set of Cu 3d and O 2p orbitals 
within the Cu04 embedded molecular cluster method and by introducing one-center (intra-center) 
Frenkel-like and two-center (inter-center) excitons. Special attention is paid to the transition matrix 
element effects both in optical and electron energy loss spectra (EELS). In the latter case we obtain 
the momentum dependence of matrix elements both for intra-center and inter-center transitions. 
We are able to give a semi-quantitative description of the optical and EELS spectra for a large 
number of OD (like CUB2O4), ID (Sr2Cu03) and 2D (like Sr2Cu02Cl2) insulating cuprates in a 
unifying manner. By comparing our analysis with the experimental data we find that the CT 
gap in insulating cuprates is determined by nearly degenerate intra-center localized excitations and 
inter-center CT excitons. The former are associated with a hole CT transition b\g — > e„(7r) from the 
b\g oc d^2_y2 ground state of dominantly Cu 3d^2_y2 character to a purely oxygen dominantly O 2p^ 
state localized on one CUO4 plaquette, whereas the latter correspond to a big — > big CT transition 
between neighboring plaquettes with the Zhang-Rice (ZR) singlet being the final two-hole state. 
The corresponding EELS intensity is found to be strongly k dependent: even for isolated exciton it 
sharply decreases by approaching the Brillouin zone (BZ) boundary. It is shown that the interaction 
of two-center excitons can result in destructive interference effects with an intensity compensation 
point. 



I. INTRODUCTION 

The nature of the electron-hole excitations in par- 
ent quasi-2D cuprates such as La2Cu04, Sr2Cu02Cl2, 
YBa2Cu306, and their ID counterparts like Sr2Cu03, 
Sr2Cu02, Li2Cu02 represents an important challenging 
issue both for the high-Tc problem and, more generally, 
for strongly correlated oxides. The mechanism of low- 
energy electron-hole excitations as well as those with 
higher energy and high intensity is still unclear. 

It is now widely believed that the most intensive low- 
energy electronic excitations in quasi-2D insulating cop- 
per oxides correspond to the transfer of electrons from O 
to Cu in the Cu02 layer, hence these materials are charge 
transfer (CT) insulators.^ Moreover, sometimes the in- 
tensive band near 2 3 eV is naively supposed to be the 
only representative of the O 2p-Cu 3d charge transfer, 
and higher-lying structures are assigned to transitions 
relevant to Cu 4s, La Sd/Af, Sr 5d, or Ca Ad depend- 
ing on the actual chemical formula of the corresponding 
cuprate. 

Despite the giant (^ 10*) number of experimental 
and theoretical papers on cuprates we deal actually 
with a lack of detailed studies of electron-hole excita- 
tions both in parent and high-Tc cuprates. The most 
part of optical information^ was obtained to date by re- 
flectivity measurements followed by a Kramers-Kronig 
transformation usually accompanied by a number of un- 
avoidable uncertainties as regards the peak positions 



and intensities of weak spectral features. Nevertheless, 
applying different optical measurements including Ra- 
man spectroscopy and nonlinear susceptibility, one ob- 
tains important information, in particular, as regards 
the subtleties of the optical gap states in such insulat- 
ing cuprates as CuO,2i4La2Cu04(SiSiLaR2Cu04 (R=La, 
Nd, Eu, Gd)ri^ YBaaCusOgjia-ii SraCuOaClafi^iiiili^ 
and SrzCuOgiSJiia 

Unfortunately, any conventional optical technique is 
momentum restricted, and yields an indirect information 
on the dielectric function e(w, k) only at the F-point given 
by k — !■ 0. In other words, such a technique cannot probe 
the dynamics of electron-hole excitations. 

In this connection, we would like to emphasize the de- 
cisive role of direct EELS measurements in the obser- 
vation and analysis of electron-hole excitations as com- 
pared with conventional indirect optical data. The pos- 
sibility to yield the polarization and momentum de- 
pendent loss function Im{—l/e{uj,'k)) makes EELS a 
powerful tool to examine subtle details of the energy 
spectrum and dynamics of electron-hole excitations in 
cuprates. ^^■^^'^"'^ In the limit k ^ the selection rules 
are the same as in optics, i.e. only dipole transitions 
are allowed. For finite k nondipole transitions are de- 
tected as well. In contrast to reflectivity measurements, 
the EELS in transmission is a surface insensitive tech- 
nique. Recently, momentum resolved resonant inelastic 
X-ray scattering (RIXS) measurements^SiS^ performed 
for Sr2Cu02Cl2 and Ca2Cu02Cl2 have demonstrated the 
feasibility of RIXS to study electron-hole excitations in 
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insulating cuprates. 

The theoretical analysis of optical excitations beyond 
simple band models has been mostly performed in the 
frame of conventional Hubbard-like models, where one 
usually considers only O 2p„ orbitals, and neglects O 
2p7r and O 2pz orbitals at all. In the parameter regime 
appropriate for undoped cuprates, one implies that the 
one hole per unit cell is mainly localized on the Cu site. 
Application of the current operator onto this state pro- 
vides a charge transfer from the Cu to the O site. Thus, 
one would expect a dominant absorption feature at the 
energy of the charge transfer gap. In terms of the Hub- 
bard model, this is a CT transition from the nonbonding 
oxygen to the upper Hubbard band. It was emphasized 
that instead of a transition to the upper Hubbard band, 
it is more reasonable to speak of a transition to the cor- 
related Zhang-Rice2i type statesiS^ 

One of the central issues in the analysis of electron- 
hole excitations is whether low-lying states are comprised 
of free charge carriers or excitons. A conventional ap- 
proach implies that if the Coulomb interaction is effec- 
tively screened and weak, then the electrons and holes are 
only weakly bound and move essentially independently as 
free charge-carriers. However, if the Coulomb interaction 
between electron and hole is strong, excitons are formed, 
i.e. bound particle-hole pairs with strong correlation of 
their mutual motion. In practice, many authors consider 
excitons to consist of real-space configurations with elec- 
trons and holes occupying nearest neighbor sites, while 
the electrons and holes are separated from each other in 
the conduction-band statesiS& 

The electron-hole excitations near the CT gap in insu- 
lating Cu02 planes were theoretically examined in Ref. 
Wh in terms of a six-band model and using the cell- 
perturbation method. The excitons were approximately 
treated as electrons and holes moving freely in their re- 
spective quasi-particle bands, except in the nearest neigh- 
borhood, where they feel an attraction and an on-site in- 
teraction determined by the eigenstate of the cell Hamil- 
tonian, which leads to a rather small exciton dispersion. 

However, a small exciton dispersion is in contradic- 
tion to the results of the EELS measurements for the 
2D model cuprate Sr2Cu02Cl2>i^ These measurements 
stimulated the elaboration of a simple CT exciton model 
by Zhang and Ng>i2i2^ They used a local model to study 
the formation and the structure of the low-energy charge- 
transfer excitations in the insulating Cu02 plane restrict- 
ing themselves to the Cu 3da,2_j^2 and O 2p„ orbitals. 
The elementary excitation is a bound exciton of spin sin- 
glet, consisting of a Cu^ and a neighboring Zhang-Rice 
singlet-like excitation of Cu-0 holes with a rather large 
dispersion. They considered four eigenmodes of excitons 
with different symmetry Aig, Big, E^, or 5, D, -Pi, 2- The 
experimental peak analysis in Sr2Cu02Cl2 revealed a sur- 
prisingly large energy dispersion « 1.5 eV along [110] 
which was interpreted in terms of a small exciton mov- 
ing through the lattice freely without disturbing the anti- 
ferromagnetic spin background, in contrast to the single 



hole motion>i2i2& So, it seems that the situation in an- 
tiferromagnetic cuprates differs substantially from that 
in usual semiconductors or in other bandlike insulators 
where, as a rule, the effective mass of the electron-hole 
pair is larger than that of an unbound electron and hole. 

Other theoretical works in the frame of conventional 
Hubbard-like models, without the O ip-j^ and 2pz orbitals, 
are those by Hanamura et. al.^^ and Kuzian et. alJ^^ The 
excitonic cluster mode?^ has been elaborated to describe 
the CT transitions in insulating cuprates in frame of the 
conventional three-band Hubbard Hamiltonian. The au- 
thors made an attempt to treat both the bound and un- 
bound states of the CT electron-hole pair on the same 
footing. To this end they performed perturbation calcu- 
lation for the excited states of Eu symmetry generated 
by the oxygen O 2pa^ hole separated from the Cu ?>dx^~y^ 
electron by distances of 1, 2, 3, 4, or 5 nearest neighbors. 
To evaluate the dipole transition matrix element they 
took into account only the transition between the dx2_y2 
and the O 2pa- orbital of e,j symmetry. The frequency 
and wave number dependence of the dielectric function 
e(k,a;) and its inverse e~^(k, has also been studied 
in Ref. [s^- The authors show that the problem, in gen- 
eral, cannot be reduced to a calculation within the single 
band Hubbard model, which takes into account only a 
restricted number of electronic states near the Fermi en- 
ergy. The contribution of the rest of the system to the 
longitudinal response is essential for the whole frequency 
range. With the use of the spectral representation of the 
two-particle Green's function they show that the problem 
may be divided into two parts: into the contributions of 
the weakly correlated and the Hubbard subsystems. For 
the latter an approach is proposed that starts from the 
correlated paramagnetic ground state with strong anti- 
ferromagnetic fluctuations. The method is applied to the 
multiband Hubbard (Emery) model that describes lay- 
ered cuprates. 

Already a shorthand inspection of the original ex- 
perimental data,^^ and especially of the high-resolution 
EELS spectraSSiSi for Sr2Cu02Cl2 point to some es- 
sential shortcomings of the Zhang-Ng (ZN) model. 
The high-resolution EELS spectra^SiSi for insulating 
Sr2Cu02Cl2 with a momentum range along [110] broader 
than in Ref. il£ did not conflrm two principal predictions 
of the ZN-model: no spots of the predicted D-mode ex- 
citon were found in a broad energy range below 2 eV; a 
clear decrease contrary to the predicted increase of EELS 
intensity along [100] was observed with momentum in- 
crease from k = 0.5 to k = 0.7. Moreover, the high- 
resolution measurementsSSiSi unambiguously detected a 
" multi-excitonic" nature of EELS spectra in the 2-^8 eV 
energy range with a clear manifestation of several dis- 
persionless narrow bands with energies near 2.0,4.2,5.4, 
and 7.2 eV. Thus, we conclude that the simple ZN model 
should be revisited to consistently explain the large body 
of new experimental data. 

Especially important is the proper inclusion of the O 
2p7r orbitals into the theory of electron- hole excitations. 
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One of the first microscopic models for the low-energy 
CT transition in the insulating copper oxide CuO peaked 
near 1.7 eV was given in Ref. This transition was 
assigned to the dipole-allowed CT transition big — > e„ 
from the Cu 3d-0 2p hybrid big oc (1^2 _y2 hole ground 
state to the purely oxygen orbital doublet state with 
predominant O 2pT^ weight. To the best of our knowledge, 
it was the first indication of an optical manifestation of 
low-lying nonbonding O 2pt^ states. 

In the present work we generalize the ZN-theory into 
a simple, physically clear, cluster theory for the CT ex- 
citons that on the one hand catches the essential physics 
and on the other hand provides a semi-quantitative de- 
scription of the optical and EELS spectra. A prelimi- 
nary analysis^i*^ has shown that a concise interpreta- 
tion of the experimental EELS spectra may be obtained 
by developing further the ZN-model with respect to two 
points. First of all, we include the complete set of Cu id 
and O 2p orbitals of the plaquette in terms of the em- 
bedded cluster method. And second, we argue that the 
elementary 0-Cu CT processiSiS^ generates two types of 
excitons: one-center excitons localized on one CUO4 pla- 
quette and two-center excitons extending over two CUO4 
plaquettes. We base our argumentation mainly on the k 
dependence of the corresponding EELS intensity as re- 
vealed by a detailed analysis of matrix elements of one- 
and two-center excitons. We also investigate the exciton 
dynamics. Finally we arrive at a semi-quantitative un- 
derstanding of the electronic states that dominate the op- 
tical and EELS response in a rather broad energy range. 
In such a way we explain the main features of optical 
and angle-resolved EELS spectra for a large number of 
different cuprates. They reach from the OD compound 
CUB2O4, over the ID cuprate Sr2Cu04 up to the 2D 
model compound Sr2Cu02Cl2 and the parent cuprates 
of high-Tc superconductors. For instance, we present an 
alternative explanation of the angle-resolved EELS spec- 
tra of Sr2Cu02Cl2 that is free from many shortcomings 
of the ZN-model. We should emphasize that our theory, 
like also those of Zhang and Ng cannot provide a proof 
whether the low-lying particle-hole pairs are bound or 
not. It has to be expected that at least the higher lying 
particle-hole excitations on a finite cluster do in reality 
not exist as quasi-particle excitations with infinite life 
time. Nevertheless, they should be visible in the spectra 
as resonances, having possibly a rather large broaden- 
ing. In that sense, we will denote all particle-hole exci- 
tations as excitons in the present paper. The question 
of its life time has to be answered independently. There- 
fore, we cannot contribute for example to the ongoing 
discussion about the stability of long-lived excitons in 
ID cuprates.Si24iS 

The rest of the paper is organized as follows. In Sec. 
II we shortly address the electronic structure and the en- 
ergy spectrum of a CUO4 cluster. The one-center small 
excitons with appropriate transition matrix elements are 
considered in Sec. III. The two-center small excitons and 
accompanying issues of correlations, final state effects. 



EELS transition matrix elements, even and odd excitons 
are addressed in Sec. IV. In Sec. V we address the prob- 
lem of the energy and intensity dispersion of two-center 
excitons. Sec. VI is devoted to a model analysis of the 
experimental EELS spectra in OD, ID and 2D insulating 
cuprates. 

II. ELECTRONIC STRUCTURE OF 
COPPER-OXYGEN CLUSTERS 

The electronic states in strongly correlated cuprates 
manifest both significant correlations and dispersional 
features. The dilemma posed by such a combination is 
the overwhelming number of configurations which must 
be considered in treating strong correlations in a truly 
bulk system. One strategy to deal with this dilemma is 
to restrict oneself to small clusters, creating model Hamil- 
tonians whose spectra may reasonably well represent the 
energy and dispersion of the important excitations of the 
full problem. Naturally, such an approach has a number 
of principal shortcomings, including the boundary condi- 
tions, the breaking of local symmetry of boundary atoms, 
and so on. 

As an efficient approach to describe excitonic states, es- 
pecially with small effective electron-hole separation we 
propose here the embedded molecular cluster method. In 
the present context we use one or two neighboring CUO4 
clusters embedded into the insulating cuprate. This 
method provides both, a clear physical picture of the 
complex electronic structure and the energy spectrum, 
as well as the possibility of quantitative modelling. Es- 
kes et al^, as well as Ghijsen et al^ have shown that in 
a certain sense the cluster calculations might provide a 
better description of the overall electronic structure of in- 
sulating copper oxides than band-structure calculations. 
They allow to take better into account correlation effects. 

Beginning from 5 Cu 3d and 12 O 2p atomic orbitals 
for Cu04 cluster with D^h symmetry, it is easy to form 
17 symmetrized ai^, 02^, 6ig, 629, (gerade=even) and 
a2ti, ^2m, e„(tT), eu(7r) (ungerade=odd) orbitals. The even 
Cu 3d aig{Mz2),big{M^2_y2),b2g{idxy),eg{3dooz,'idyz) 
orbitals hybridize, due to strong Cu 3d-0 2p covalency, 
with even O 2j5-orbitals of the same symmetry, thus 
forming appropriate bonding 7** and antibonding 7° 
states. Among the odd orbitals only eu(iT) and e„(7r) 
hybridize due to nearest neighbor pp overlap and trans- 
fer thus forming appropriate bonding and antibonding 

purely oxygen states. The purely oxygen a2g, a2u, b2u 
orbitals are nonbonding. All "planar" O 2p orbitals (see 
Fig. in accordance with the orientation of lobes could 
be classified as a {aig,big,eu{cr)) or tt (a2g, 623, e„(7r)) 
orbitals, respectively. 

Bonding and antibonding molecular orbitals in hole 
representation can be presented as 

I7'') = cosa^|7(3d)) -I- sina^|7(2p)), 
I7") = sina-y|7(3d)) — cosa-y|7(2p)). 
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FIG. 1: Electron (hole) density distribution for planar copper 
and oxygen molecular orbitals of one Cu04 plaquette. 



For example, 

l^ig) = cosafciJ6ig(3d)> +sinah,J6ig(2p)), 

\btg) ==sinabiJ6ig(3d)) -cosaf,iJ6ig(2p)), (1) 
where big{3d) — ?>d^2_y2\ 

|e^) = cosofe |e„(7r)) + sinag |e„(CT)), 



|e°) = sinofe |e„(7r)) - cosag |en(cr)) 



(2) 



equally for both types (a;, y) of such orbitals. 

To explain the electronic structure of the embed- 
ded molecular cluster in a broad energy range we need 
the complete information about the bare parameters 
of the effective Hamiltonian. The numerical values of 
energy parameters for holes in insulating cuprates like 
La2Cu04 or Sr2Cu02Cl2 which were used by different 
authors^SiSMSiiii^ can be summarized as follows (in eV) : 

ep — ~ 2 

Ud^ 6^10 eV; C/j,w4h-6; Vpd w 0.5 1.5; 



E(eV) 




bai-cCu3d ^g-ab7^ 



- 



\ 

Cu3d-CF effect 3d-2p 

hybridization 



2p-2p 
,' hybridization 



02p-CF effect 



- 



FIG. 2: Model single-hole energy spectra for a CUO4 plaque- 
tte with parameters relevant for Sr2Cu02Cl2 and a number 
of other insulating cuprates. 



ipd ~ 1-0 . 1.5, l^ppcrl ~ 0.8, ippTT ~ r)l^ppcr| ^ 0.4, 



where 



1 _ 1 

T / ^ ^d,i 0, Cp lyX^pfy ~^ ^'PTT ^" ^pz) 



are the centers of "gravity" for the Cu 3d and O 2p man- 
ifolds, respectively. As will be discussed more in detail 
below (see Sec. VI.E.l) such a parameter choice is con- 
firmed by the corresponding ARPES datat^^iM giving im- 
portant information concerning the nonbonding oxygen 
states. Fig. |21presents a possible single-hole energy spec- 
trum for a Cu04 plaquette embedded into an insulating 
cuprate like Sr2Cu02Cl2 calculated with the parameters 
above. For illustration we show also a step-by-step for- 
mation of the cluster energy levels from the bare Cu 3d 
and O 2p levels with the successive inclusion of crystalline 
field (CF) effects, O 2p-0 2p, and Cu 3d-0 2p covalency. 



III. ONE-CENTER FRENKEL-LIKE EXCITONS 

Small charge transfer excitons, or excited states arising 
from the configuration in which an electron is transferred 
from a negative ion to a nearest neighbor positive ion 
have been considered many years ago to be the origin of 
excitonic-like peaks which show up as a low-energy struc- 
ture of the fundamental absorption band, for instance in 
alkali halide crystals. Namely this idea was exploited by 
Zhang and Ng in their simple model theory of small CT 
excitons in insulating cupratesii2i2& However, in terms 
of the Cu04 cluster model the local two-atomic Cu 3d- 
O 2p charge transfer generates a number of both intra- 
cluster and inter-cluster transitions (see Fig. The 
corresponding electron-hole excitations form small one- 
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FIG. 3; (a) Simplified illustration of the one-center exciton 
and (b) the two-center exciton of hig\pd, i.e. ZR-singlet type, 
(c) Classification of the important two-center excitons to- 
gether with an energy scheme. 



and two-center excitons, respectively. The main crite- 
ria are the binding energy and the average electron-hole 
separation. Namely the latter provides the localization 
and stability of excitations due to their small overlap 
for nearest neighbors (small tunneling amplitude), and 
finally leads to a narrow exciton band. Below we ad- 
dress the most probable candidate states for small one- 
and two-center CT excitons. It should be noted that the 
main concept which we apply has in many aspects a close 
similarity with those used for CT excitons of linear-chain 
organic 7r-conjugated polymersi2Sii& Some symmetry as- 
pects of the CT excitons in cuprates were addressed by 
Cherepanov et alm^ 



A. Classification 

Among the numerous one-center electron-hole CT ex- 
citations the first candidates for dipole-active excitons 
are CT transitions b\g e"''' from the ground state 
b\g to the purely oxygen doublet O 2p7r-0 2pcr hybrid 
states e'^'^, which are allowed in "in-plane" polarization 
E _L C4, and the b\g — > 62M transition to purely oxygen 
O 2pz-\ike state, which is allowed in " out-of-plane" po- 
larization E II C4. All these one-center excitons may be 
rather simply represented as a hole rotating on the four 
nearest oxygens around an electron predominantly local- 
ized in the Cu 3dx2-y2 state with a minimal electron-hole 
separation Reh « RcuO ~ 2A (see Fig.lJJ. 

The two excitons h\g — > and h\g e° differ by the 
oxygen hole density distribution: for the former this has 
a predominantly O 2pTr character, while for the latter it 
has a O 2p„ one. It should be noted that the excitonic 
doublet state may be considered as two current states 



eti(±l), or two currentless states e„(a;,?/) with quenched 
orbital motion. The former are excited by circular polar- 
ized light, while the latter by linear polarized one. The 
double degeneracy for the b\g e\ excitons results in 
anomalously strong electron-lattice coupling with all fea- 
tures typical for the Jahn- Teller effect^ in particular, 
high probability to form a self-trapped state. In any case, 
these excitons should be considered as first candidates for 
resonant phonon Raman scattering activity. One should 
specifically emphasize that along with resonance excita- 
tion of the allowed Raman modes the Jahn- Teller exci- 
tons could generate forbidden phonon modes. 



Among dipole inactive excitons one should note h\g — > 
b2g and, especially, b\g a2g with purely oxygen O 2p7r 
holes. The latter incorporates a purely oxygen a2g hole, 
which has in frames of our model the minimal energy 
among all the purely oxygen states. In accordance with 
the model energy spectrum (see Ref. lsif) its energy should 
be of the order of 1.5-;-1.8 eV, or in other words appears to 
be lower than the optical gap. This circumstance draws 
specific attention to this exciton, despite the b\g —> a2g 
transition is dipole forbidden. Due to the purely oxygen 
nature of the hole, the b\g — > a2g exciton resembles the 
dipole active — *■ e„ exciton, however, contrary to the 
latter this is a purely O 2p.^ hole. 



There are several examples of a qualitative assignement 
of low-lying optical spectral features to the non-bonding 
oxygen orbitals, appearing in CuO near 1.7 eV (see Ref. 
m and in SrzCuOsCls near 2.5 eV (see Ref. El). Also, 
it was shown that the dipole-allowed one-center electron- 
hole excitation b\g — > e„(7r) is visible in the low lying part 
of the EELS spectrum of Sr2Cu02Cl2 However, to the 
best of our knowledge, there is at present no unambigu- 
ous identification of this one-center excitation (and also 
of big b2u) as well-separated entities in 2D insulating 
cuprates. That is different for ID cuprates where these 
one-center excitations built up of non-bonding oxygen- 
orbitals could be separated in EELS by choosing a mo- 
mentum transfer perpendicular to the chain directioniS 
The conventional one- and two-band Hubbard-like mod- 
els consider only the Cu 3dj.2_y2 and O 2p^ orbitals, and 
do not take into account the nonbonding O 2p7r and O 
2pz orbitals. Partly, that may be justified by the specific 
properties of the big {dr^2_y2) ground state resulting in a 
predominant contribution of O 2prj orbitals both to the 
Cu 3d - O 2p bonding and to the different charge transfer 
transitions. Optical absorption and EELS activity asso- 
ciated with the nonbonding O 2p-^ and O 2pz orbitals 
appear to be rather weak if any. 
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B. Transition matrix elements in optics and EELS 

1. Intensities of the electric dipole CT transitions for 
one-center small excitons 

Electric dipole CT transitions for the hole localized in 
the Cu04 cluster are allowed from the even ground state 
b\g to the bonding and antibonding purely oxygen odd 
e"''' states for the E _L C4 polarization, and to purely 
oxygen odd 62M state for the E || C4 one (see Fig. 
Dipole transition matrix elements for the b\g ground state 
depend essentially on the fact whether the final state is of 
O or O 2pcr character. For the former case we have 
only a "nonlocal overlap contribution" with two-center 
integrals 



cos {d^2^y2\qx\eux{T^)) 



1V2 



cos ttbig 



{blg\qz\b2u) = qV2cosab,^ J <j)d^2_y2i^)* z<Ppzir - R)dr, 



where R || y, and 



(3) 



that results in equal intensities for b\g eu{TT) and 
b\g b2u transitions. Naturally, that in a more common 
case we have orthogonalized Wannier functions instead of 
simple atomic ones, and the calculation of the nonlocal 
overlap matrix elements is not so straightforward. For 
the b\g e„((T) dipole transitions to purely O 2pa- state 
we have both nonlocal and local ( or " covalent" ) overlap 
contributions due to nonzero O 2p^ density in the initial 



state 



(b'lglqxleuxicr)) ^ -(6jgk2/|e„a(cr)) 



cos ttfcig 



:ir)*y(jjpy(r -R)dr 



qRcuc 



■ sm afeig 



where R || y. It should be noted that namely the lo- 
cal term is the only contribution usually addressed for 
transition matrix elements. 



The bonding and antibonding e"''' states incorporate 
both O 2p„ and O orbitals. Therefore, the respective 
dipole transition matrix elements acquire a rather com- 
plex structure. However, for the "covalent", and proba- 
bly, the leading contribution, we obtain a simple relation 



I{b 



= tan ttp 



(4) 



In other words, the relative intensity of the two main 
elcctric-dipole transitions is mainly determined by the 
magnitude of pp- hybridization and does not exceed 0.1 
for the typical values of the parameters. 

In some cases, it is of practical interest to provide in- 
formation concerning the allowed dipole transitions be- 
tween excited states. In our case, these are, first of all, 
numerous transitions from or to e„ states. For instance, 
the a2g — ^ dipole allowed transition has a very large 
intensity due to a large local contribution to the matrix 
element: 



{a2g\qx\eux{'^)) = 



qR. 



CuO 



^/2 



It should be noted that the rather simple model of a 
single CUO4 cluster allows already to predict both the 
energies and the relative intensities for such transitions. 



2. Intensities of the EELS transitions for one-center small 
excitons 

The intensity of an exciton is determined by the imag- 
inary part of e(k, lo) around a pole, and is given by 



(5) 



where "^exc is the exciton wave functioniS^ Within the 
long-wavelength approximation {k — > 0) the EELS selec- 
tion rules are the same as in optics if we address the k 
direction to be that of the electric field. For the one- 
center small 7-excitons it is rather easy to obtain the 
angular k-dependence of the transition matrix elements 
in the small k-approximation: 

(vE'e.c|e''"'|*Gs) (X (7|e^'''-|6?,) (x fc^lL^(k), (6) 

where YL^{)<i) is a linear superposition of spherical har- 
monics YLM(k), which forms a basis of the irreducible 
representation 7 of the D^h point group. We make use 
of the familiar expansion 



(7) 



L=0 M= 



restricting ourself to the first nonzero contribution with 
minimal L. Thus, for different 7 one obtains: 

(e„|e''"'|6i ) cx fcYie^(k) cx kx^y', (dip. approx.) 
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(aig|e''"'|&?g> « fc^F26i<,(k) cx (fc^ - fc^); (quadr. approx.) 
(a2g|e*'''"|65g) oc fc^r2fc2g(k) oc k^ky-, (quadr. approx.) 



IV. TWO-CENTER EXCITONS 
A. Electronic structure and classification 



{big\e^^'^\big) oc k'^YiaigO^) « cosAcj); (oct. approx.) 

(^2ff|e''"'|6ig) cx k^Y4a2gOi-) oc sin 40; (oct. approx.) . 

Here, (p is the azirauthal angle of orientation for the k 
vector. In general, among dipole inactive excitons for 
corner-sharing CUO4 plaquette systems the aig and a2g 
ones are "visible" only along [100] and [110] directions, 
respectively; the b^g exciton is visible both along [100] 
and [110] directions, while the b2g one does not mani- 
fest itself both along [100] and [110] directions. Taking 
into account only the conventional local overlap ("cova- 
lent" ) contribution to the transition matrix element one 
may obtain model expressions for {'^\e^^'^\b\g) suitable 
throughout the BZ. Thus, one obtains for the O 2p local 
overlap contribution: 



{a2g\e^''\b\g) 



'is 



(b 



a I „ikr i l6 



19 1 



1^* 



19/ 



(etxle 



smabj^ cosQfaig 



■ sma^jg cosabig 



19/ 



cos ■ 



0; 



cos 



ky (X 



COS ■ 



cos 



(8) 



V2' 



kxd 



{eiy\e''"'Kg) 



[e„x\e 



Ky\e'nb\g) 



% , . . kyCt 

1= sm atj^g sm ae sm — ^ ; 



— =sinah, sinaeCos— 
V2 2 



1 , . kyOi 

-^=sma6, cosQfeSm— 
V2 " 2 



For the dipole- allowed big e„ transitions the EELS 
intensity has to decrease in going from the F-point to the 
BZ boundary with 

/(7r/a,Va,0) ^ /(7r/a,0,0) _ 4 ^ 
7(0,0,0) 7(0,0,0) 7r2 ~ ■ ■ 

One should note that with the inclusion of the copper 
contribution the final expression for the matrix element 
(&?g|e*'""|65g) should be modified 

(&?,|e*"-|&L) 



1 

--smabig cosckfcig 



cos ■ 



kyU 

cos — 2 

2 



(9) 



Interestingly, that in the frames of the local model for 
transition matrix elements the a2g and one-center ex- 
citons appear to be invisible throughout the BZ body. 



Inter-center charge transfer transitions between two 
CUO4 plaquettes centered at neighboring sites A and B 
define two-center excitons in the molecular cluster CU2O7. 
These two-center excitons may be considered as quanta 
of the disproportionation reaction 



CuOr + CuOr ^ CuOr + CuOr (lO) 



with the creation of electron CUO4 and hole CUO4 



cen- 



ters. The former corresponds to completely filled Cu 3d 
and O 2p shells, or the vacuum state for holes, while the 
latter may be found in different two-hole states. These 
two- hole states can be classified by 7172 = F according 
to the symmetries % of each of the holes. For instance, 
if we restrict ourselves to charge transfer processes gov- 
erned by the strongest a bond, one has to distinguish 
three different channels bfg, bigaig and 5ige„. Only two 
of them (big and bigCu) have non zero optical matrix 
elements, but all arc important for EELS. Each symme- 
try different channel contains several wave functions with 
different Cu and O character (see Section IV. A. 1). 

To construct ground state and charge transfer excited 
states we use an approach similar to the well known 
Heitler-London scheme. In this approach the two hole 
states for the CU2O7 molecular cluster are formed from 
different one-center states. We deal with charge states of 
the CUO4 plaquette with zero, one, and two holes, respec- 
tively, denoted as Then, the ground state "^gs 
of the cluster CU2O7 is predominantly given by the sym- 
metrized product of *&^^(feig)'I'B^(6ig) coupled to higher 
states of the same symmetry. Accordingly, the excited 
(exciton) states are built by product states 



'/'f5(r) = <i>i°^$g'(r) 



(11) 



with subsequent (anti)symmetrization (see Section 
IV.A.2). 



Two-hole configurations for the CuO^ center 



To classify all two-hole configurations of one CUO4 
center we use their symmetry F = 7172 and take into ac- 
count the Coulomb interaction. The latter is of particu- 
lar interest both for energetics and the electron structure 
of two-center excitons. According to symmetry we have 
to distinguish even F configurations (like bfg or bigGig) 
and odd F configurations (as bigCy). Alternatively to the 
notation 7172 we will also use the representations of the 
7)4/1 symmetry of the CUO4 cluster (see Ref. 36) with the 
following correspondence: b\g — ^Aig and bigCu — ^Eu- 
Each channel of two-center excitons, characterized by a 
7172 configuration, consists of several wave functions dif- 
fering by their partial Cu 3c? and O 2p density distribu- 
tion. These excitons will be denoted by 7172; dd, 7172; pd. 



8 



and 7172; pp, respectively. In the sense of quantum chem- 
istry the mixing between several wave functions of one 
channel can be understood as a configuration interaction 
effect. 

Let us present the simple example of the calcula- 
tion of the two-hole spectrum in the Zhang-Rice (ZR)- 
singlet sector^^ (i.e. in the 6fg-channel). As usually, we 
assume that the ZR-singlet represents the lowest spin- 
singlet state formed by the interaction of three "purely 
ionic" two-hole configurations \d?), \pd), and Here, 
\d) — \dj;2_y2) and \p) — \pbig) are the non- hybridized 
Cu 3dx2-y2 and O 2pa- orbitals, respectively, with bare 
energies Cd and Cp. Then, the matrix of the full effective 
Hamiltonian within the bare basis set has a rather simple 
form 

/2ed + Ud t \ 

H=\ t ed + e.p + Vpd t , (12) 

V t 2ep + u;/ 

where the effective Coulomb parameter for purely oxygen 
configuration incorporates both the intra-atomic param- 
eter Up and the oxygen-oxygen coupling to the first and 
second nearest neighbors, respectively 




and the following condition holds: 

Ud > Up > Vpd 

For reasonable values of parameters (in eV): Ud — 8.5, 
Up = 4.0, Vpd = 1.2, ed = 0, ep = 3.0, t = tpd = 1.3 (see 
Ref. li^ we obtain for the ZR-singlet energy Ezr = 3.6, 
and its wave function 

= \big;pd) = -0.25|d2)+0.95bd)-0.19|p2) . (13) 

It reflects the well-known result that the ZR-singlet rep- 
resents a two-hole configuration with one predominantly 
Cu 3d and one predominantly O 2p hole. The two ex- 
cited states with energy Ezr + 5.2 and Ezr + 6.7 eV are 
described by the wave functions 

|$(')) = \blg; dd) - 0.95|d^) + 0.21|pd) - 0.22^2), (14) 

1*3"^) = \blg\PP) = -0.17|d2)-0.24bd)-0.96|/), (15) 

respectively. Given the ZR-singlet energy one may cal- 
culate the transfer energy from h\g to the neighbor- 
ing blgipd state (the ZR-singlet state), or in short the 
b\g big transfer energy: 

Act - Ezr - 2Eb,^ = (3.6 + 0.5)cV = 4.1eV, 

where the stabilization energy for the bonding b\g state is 
simply calculated from matrix l(T^ at Ud — U* — Vpd — 
0. One should note that the minimal b\g big transfer 
energy relatively weakly depends on the value of the pd 



transfer integral; we obtain the same value 4.2 eV at any 
value of t between 1.0 1.5 eV. This energy depends 
mainly on the values of (e^ — Cd) and Vpd- 

Similar charge transfer energies are obtained in various 
model approaches. It should be emphasized that this 
quantity plays a particular role as the minimal charge 
transfer energy which specifies the charge transfer gap. 
In the general case it is defined 

Act — En+i + i?jv-i — 2En, 

as the energy required to remove a hole from one region 
of the crystal and add it to another region beyond the 
range of excitonic correlations. The exact diagonalization 
studies for a series of clusters with different siz o^^i^° show 
that Act strongly diminishes with cluster size from « 
4eV for small clusters to 2.5eV as extrapolated value 
for large clusters. 

2. Ground state and charge transfer excited states of the 
Cu2 O7 cluster 

As already noted, the states 'I'J^'jj'^' of the CUO4 pla- 
quette centered at A or B will be used to construct the 
CU2O7 wave functions in a Heitler-London like scheme. 
Moreover, for two-center molecular clusters like CU2O7 
with D2h point symmetry one has two types of excita- 
tions: 0fs(r) = <i>o,<i>|(r), and </>^'=s(r) = $^(r)$^, 
which differ only by the permutation of electron and hole. 
These functions will interact due to the resonance reac- 
tion 

CuOl" + CuO^" ^ CuO^" + CuOl" , (16) 

giving rise to the splitting of the two bare excitations 
discussed in more detail below. 

The spin-singlet even ground state for the CU2O7 sys- 
tem can be written as follows 

*GS cos a 0^, , +sina Var0|s(r), (17) 
where (f)^^ ,5 is the symmetrized product of the one- 

Ig Ig 

center one-hole states ^A\big)^B\big)- coupled 
to the even states 

0|<,(r) = ^(0f5(r) ± 4s(r)) (i8) 

corresponding to symmetrized (antisymmetrized) bare 
states (j)ll^c;{T) for even (odd) F. Correspondingly, the 
odd and even excited CT states are given by 

r 

*ls = cos a V'cr (F) - CSS sin 6 , (19) 
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where the index T labels different final configurations 
of the two- hole CUO4 cluster denoted by F = 7172; dd, 
Jij2',pd, and 7172; pp, respectively. In close analogy to 
(|18|l we can write 

r^si^) = ^(^fs(r) T d^fsim (20) 
for even (odd) T. 

3. The S- and P-Uke two-center excitons 

The energies of the superposition states (|18|l and H2U|) 
are given by Eq ± |Tas|, where Eq is the energy of the 
bare state (F) and Tab is an effective resonance two- 
particle transfer integral of the resonance reaction 
The even (odd) states {(t^Es) ^^^^ denoted as 
(Biu) states corresponding to the D2h point symme- 
try of the CU2O7 cluster. Alternatively, they correspond 
to S- (for Ag) or P-like (for Si„) two-center excitons. 
Namely these S- or P-like excitons will draw our main 
attention below. For brevity, we use the simple g, u labels 
instead of the Ag, _Bi„ notations. 

The different coefficients in (|17f) and (|19|l are coupled 
by normalization and orthogonality conditions like 

XI l^rP = XI l^rP = 1; CES = ^a^cr. 
TV r 

Its magnitude is actually determined by appropriate hole 
transfer integrals. It should be noted, however, that our 
approach is only an approximative one. So, due to the 
small overlap between the two neighboring CUO4 plaque- 
ttes, the one-hole 7 label is likely to be a "bad quantum 
number". Nevertheless, it can be used in the frames of 
a semi-quantitative approach. In the single channel ap- 
proximation where the exciton is dominated by the trans- 
fer from h\g to one neighboring &ig7 state, we may write 



where i^^^^^ and Afji^^-y are the CT transfer integral 
and CT energy, respectively. 



Let us note that in our approach the S- and P-excitons 
are centered at the central oxygen ion of the CU2O7 clus- 
ter. That is in contrast to the simple ZN-model2& which 
implies the hole location on one of the CUO4 plaquettes. 
In other words, the ZN-model treats the electrons and 
holes in the e — h pair substantially asymmetrically, that 
does not allow to introduce the S- and P-like excitons. 

The magnitude of the effective resonance two-particle 
transfer integral Tab which determines the even-odd 
splitting is of particular interest in exciton theory. In the 
real cuprate situation the bare two-center exciton repre- 
sents a system of two neighboring electron Cu04~ and 
hole Cu04~ centers (given common oxygen). Thus, the 
resonance reaction corresponds to inter-center transfer of 
two holes, or two electrons. 

The S-exciton is dipole-forbidden, in contrast to the P- 
exciton, and corresponds to a so-called two-photon state. 
However, these two excitons have a very strong dipole- 
coupling with a large value of the S-P transition dipole 
matrix element. This points to a very important role 
played by this doublet in nonlinear optics, in particular 
in two-photon absorption and third-harmonic generation 
effectsAi& 

The large dipole matrix element between excited S- 
and P-excitons 

d=\{S\d\P)\^2eRcuCu, (22) 

is important for the explanation of the nonlinear optical 
effects in Sr2Cu03iiSii& The magnitude of this matrix el- 
ement yields a reliable estimate for the effective "length" 
of the two-center CT exciton. 



B. Transition matrix elements for two-center 
excitons 

1. General expressions 

In general, the expression for the EELS transition ma- 
trix element for two-center excitons has a rather compli- 
cated form 



„zkr I 



ES\ 



^-GS) 



sm a cos a 



-CES{ 



In? 6? 



rr' 



6|s(r)|e^''''l0|;s(r )) 



(n5|e^'^""l*G5)-sina^5>y(</.^5(F)|e^'^"-|*ls(r')). (23) 
rr' 

I 



It should be emphasized that it is the second term in important for the local contribution to the transition ma- 
the ground state wave function (fTTjl which is particularly 
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trix elements. The different transition matrix elements 
for the "two-plaquette" CU2O7 system can be easily re- 
duced to " one-plaquette" two- hole matrix elements. One 
has to use (|llll and H18I20|I which give 0^^(r) as super- 

I 



for r and T' of the same parity. Otherwise, the right 
hand sides in H24|l should be interchanged. At the F- 
point k = all these matrix elements are simply reduced 
to appropriate overlap integrals. The GS g-state tran- 
sitions at this point are symmetry forbidden, while the 
GS — > u-state transitions are allowed in the dipole ap- 
proximation (k — > 0). It is of particular importance 
to note that in accordance with the relations H24|l the 
nonzero dipole response at the F-point is derived only 
from either diagonal " one-plaquette" two- hole matrix el- 
ements, or from non-diagonal ones with different parity 



positions of the " one-plaquette" two-hole wave functions 

(2) 

$^^(F) centered at A and B. Then one finds for the 
X-axis oriented CU2O7 molecule 



(24) 

I 

of the left and right hand side's functions. 



2. Some basic transition matrix elements for CT governed 
by a bond 

Below we list some nonzero matrix elements for the 
X-axis oriented two-center S- and P-excitons generated 
by the main b\g — > aig,6ig,e„ CT transitions with the 
strongest a bonds: 



(0ls(r)|e^''-|'^l;5(r'))=cos^($(f)(r)|e^'^"-|^i'^(r')), 



One obtains for the local overlap contribution due to O 2p^ states for 7,7= aig, big 

{<t^% le^'-ne ) =PG5(02p6?,)cos^(cos^+cos^), 

(</'l5(M|e^''""l0l5(fc?.7)) = i(pi,5(02p7) +PG5(02p6?^))cos^(cos^ -l-cos^), 
{<t>'Es{^\lW''WEs^b\g^)) = \c*Es{02viyEs(P2v^) COS ^(cos ^ - cos ^), (7' + 7), 

{<\>'hs{b\gl)V''^'^%s{b\gl)) - -^(PBS(02P7) +PG5(02p&J^))sin:^(cos:^ +cos:^), 
(.<l>\s{b\,l)\e'^''WES^\Kgl )) ^ -^^^^^^ (7V 7); 

k'^'Es^b\ge^{<yW'WEs^b\g^)) = ic|,s(02pe„(a))cs5(02p7)sin^sin^; (25) 
(0^5K,e„(a))|e'''-|'/'l5(&L7)) = ^4s(02pe„(a))c£;5(02p7)sin^cos^; 
(0|;s(^le„(a))|e^''-|0|5(6?,e.(a))) = pEs^Veuio)) cos ^ cos ^; 
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and for the focal overlap contribution due to Cu 3d^2_y2 states 

Iff ig ig ig Z 

('/'l;s(^L7)|e*'"l'^|;s(^L7)> = (ps5(C^.3d7) + PG5(Czi3d&J,)) cos 
(</'^s(M|e^'""l0|;s(^?,7)> - ~^{pEs{Cu3d^)+pGs{Cu3dbl))s■m'^, 

(</'ls(^Le„(a))|e^''-|'/'l;5(&J.e«M)) =PG5(C«3d&?^)cos^, (26) 
(</'^s(&'ige„(0)|e^'"'l'^ls(^Le„(a))) = -^pGs{Cu3db\^) sin'^ , 
fe(&L7)|e^''-|<^ls(^L7')> = (0^5(&l7)|e^'^1</'|5(^l7')) -0, (7V 7)- 



Here, CGS.EsiCu^d'j), CGs,Esi02pj) and p = |cp are the 
probability amplitudes and hole densities on copper (oxy- 
gen) molecular orbitals for ground, or excited states, re- 
spectively. Making use of these expressions we can read- 
ily obtain the complete set of EELS transition matrix 
elements for dipole-allowed and -forbidden two-center ex- 
citons. Moreover, they can be used for the description of 
CT excitons in terms of purely copper 8^7, or oxygen 2p7 
one-hole states. Namely such an approach could be prob- 
ably more reasonable in the case of strongly correlated 
" one-plaquette" two-hole states. 

Formulas H23|) - (|26|) clearly show that the simple dipole 
approximation for EELS transition matrix elements like 
that one used in the ZN-modeli^iS^ cannot satisfacto- 
rily describe the k-dependence of the EELS intensities, 
and may be misleading. Interestingly, that the oxygen 
and copper contributions to the local overlap mecha- 
nism of EELS transitions have substantially different k- 
dependence. So, the O 2p contribution to both dipole- 
allowed and -forbidden transitions for 7,7 = aig,big 
turn to zero at the (tTjTt) point, contrary to the copper 
contribution. At this point the active oxygen contribu- 
tion to the loss function is associated only with dipole- 
forbidden '^gs — > 0£s(^ig^" (''')) transitions. For the 
[100] direction the oxygen contribution to the EELS in- 
tensity for the dipole-allowed transitions decreases twice 
as faster than the copper one. 

In other words, the EELS intensity at the (tt, tt) 
point reproduces in such a case the spectral distri- 
bution of the two-center excitations generated by the 
b\g — > Cu CT transition. In contrast, the EELS inten- 
sity near the F-point, where the S-exciton contribution 



I 

turns into zero, reproduces the spectral distribution of 
the two-center dipole-allowed excitations g ^ u, gen- 
erated by the b\g big CT transitions both to pre- 
dominantly oxygen {Ap^^ {02p'-f) > p^^ {Cu3dj)) and 
copper {4:p^^{02p'-f) < p^^ (CuSd'y)) states, as well as 
the b\g e„(7r), et((cr) CT transition to purely oxygen 
eti(7r), eu{<j) states. In principle, this makes it possible to 
examine the Cu 3d and O 2p partial composition of the 
excitons by means of the EELS intensity analysis, though 
it should be noted that dipole-allowed two-center excita- 
tions could interact both with each other and with the 
dipole-allowed one-center excitations b\g — > eti(7r), e„(CT). 



3. Transition matrix elements for the big -channel of 
mter-center CT excitons 

For illustration, let us consider the fe^^-channel of inter- 
center CT excitons. We have to take into account the 

(21 

admixture of all three singlet wave functions of the 
ZR-singlet sector f(ll3 |l - (|15|l ') to the ground state wave 
function ^'gs- The corresponding amplitudes in Eq. I|17|) 

will be denoted by a^^^ . So, the b^^-channel consists of 
three inter-center CT excitons with energy separations of 
AE12 « 5.2 eV and Ai?i3 « 6.7 eV, respectively. If we 
neglect for a moment the configuration interaction effects 
in the final states, the relative intensities of these three 
excitons at the F-point are given by: 
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More accurately, the intensity of the low-energy dipole- 
allowed two-center exciton with the formation of the ZR- 
singlet b\g\pd is determined by the matrix element: 



'■'ES\ 



{b'-pd)\e^^^\^GS 



kxCi 



-% sma sm ■ 



function. These amplitudes depend mainly on the effec- 
tive h\g big transfer integrals for the bare b\g hole be- 
tween two neighboring plaquettes. The parameters sug- 
gest the following relationship 



«Si > i<i > 



(l)l 



,(2) I 



(30) 



(0.06a« 



(d^le'^^ld^) {0.06a[]] - 0.24a[f + 0.04a 



+ {pdy^''\pd) (0.90a[^^ -I- 0.204^^ - 0.22a^f 

+ ip'le^^'^lp') (0.04a« + 0.04aiJ + O.lSa^)] (28) 
where the one-plaquette matrix elements are: 



(pdle^'^'^lpd) = 1 



1 k^a 
- cos — — 
2^ 2 



cos-|-), 



(p^\e^^'^\p^)^{cos'^+cos^) 



(29) 



Let us draw the attention to the different k dependence of 
the "diagonal" (the first term) and "non-diagonal" (the 
second and third terms) contributions to the overall ma- 
trix element. When one moves from the F-point to (tt, tt), 
the absolute value of the former decreases, contrary to the 
latter which grows from zero at the F-point to a non-zero 
value at the BZ boundary. The final k behavior of the 
overall transition matrix element could be rather compli- 
cated depending on the relative signs and magnitudes of 
the a[|*^ coefficients, or the admixture amplitudes of the 
different ZR-singlet-like states in the ground state wave 



In other words, the largest effective integrals will be ex- 
pected for the transfer to predominantly oxygen 6ig-like 
hole states. In such a case, we cannot exclude the ap- 
pearance of an intensity compensation point due to the 
competition of the " diagonal" and " non-diagonal" terms 
in the matrix element. In any case, we may expect a 
rather unusual k behavior of the transition matrix ele- 
ment which in turn determines the behavior of the EELS 
intensity for one of the main low-energy contributions. 

(3) 

On the other hand, the large value of al^^ leads to an 
higher intensity of the high energy CT exciton at the 
F-point in comparison to the low-energy one with an es- 
timated energy separation between the two of w 6.7 eV. 
All this illustrates the particular importance of the cor- 
relation effects for two-hole configurations being the final 
states for the inter-center CT transitions under consid- 
eration. In addition, one should note that a reasonable 
analysis of the b\g — *■ big channel points to a rather wide 
spectral range of the intensive two-center CT transitions. 

If it were possible to neglect the effects of d^- pd- and 
p^-mixing in the bfg configurations, the dipole-allowed 
channel could be divided into independent d^- pd- 



"is 



and p -contributions with the corresponding EELS inten- 
sities Id2 (X |(d2|e*r|d2)|2^/p^ cx |(pd|e''^"" |pd) | ^ , /p2 cx 
|^p2|gikr|p2^|2 ^ggg ^j^g matrix elements from ((23 )■ 

The intensities of dipole-forbidden {g — g) CT transi- 
tions in the 6f ^-channel turn into zero both at the F-point 
and at the boundary of the Brillouin zone. To illustrate 
the k-dependence of the appropriate matrix elements, we 
restrict ourselves below to the simplified model where 
only the ZR-singlet contributes to the 6jg-channel 



J 



b%s{big;pd)\e^ 



kr 



1 k a 

*Gs) = 2 {PEs{02p) - pGs{02p)) cos 



1, 

-(cos 



2 



kya 



(31) 



r 



It should be noted that for the [11] and [10] directions 
the matrix element reaches its maximum at k = '^ky^ax- 
This value depends on the actual difference in the big 
hole density distribution of ground and excited states. 



4- Transition matrix elements for the bi^aig -channel of 
inter-center CT excitons 

Taking into account the correlation and configuration 
interaction effects there appear three ^Big terms with a 



6igaig-like configuration which form the final states for 
the b^gflig-channel of the two-center CT transitions. This 
channel is only governed by the oxygen contribution, due 
to the orthogonality of the Cu 3dbig and Cu Sdaig or- 
bitals. A shorthand analysis of the k dependence for the 
different matrix elements listed above in Eqs. H25|l - H26|l 
shows that the ^JgCig-channel is "silent" both at the F- 
point and all along the (tt, tt) direction. Along the (tt, 0) 
direction we deal with an unusual behavior of the appro- 
priate P{Au) and S{Ag) excitonic modes. Going to the 
BZ boundary, the intensity of the former mode increases 
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while that of the latter one turns into zero both at F and 
at (tt, 0) with nonzero value in between. So, one may 
conclude that the b\gaig-cha.nnel could be revealed only 
along the (tt, 0) direction outside the F-point. 



5. Transition matrix elements for the b^gBu-channel of 
inter-center CT excitons 

Like the 6jgaig-channel the 6jge„-channel is only gov- 
erned by the oxygen contribution. As it was mentioned 
above, we have two types of purely oxygen orbitals: 
e„((T) and e„(7r) with a and tt directions of the O 2p 
lobes, respectively. These orbitals hybridize with each 
other owing to the O 2p-0 2p transfer. As a result we 
obtain four ^E^ terms for the 6ige„-like two- hole con- 
figurations (b\geu{'rr),b'igeu{cr),bfgeuiTT),b'[geu{cr)) which 
form the final state of the &^ge„-channel of the two-center 
CT transitions. 

Non-zero transition matrix elements in Eqs. 125|l - l|2t)|) 
are only obtained for the 6jgeu((T)-channcl. Therefore, 
the resulting magnitude of the transition matrix elements 
would be firstly determined by the eu{(y) weight in the 
bigBu configuration for the final ^Eu term. On the other 
hand, the main contribution to transition matrix ele- 
ments is determined by the purely oxygen big compo- 
nent of the final states. Both these conclusions are of 
particular importance for the relative intensities in the 
&ige„-channel. The doublet of transitions with final anti- 
bonding big state is more intensive than the similar dou- 
blet with the final bonding b\g state: its relative inten- 
sity could be estimated to be | tanoff^g ^ « 1.3. For either 
doublet the strongest intensity is predicted for transitions 
with final e„((T) state with estimated relative intensity to 
be I tanae„ P > 1. 

Different theoretical estimations (see above) point to 
the predominantly 6ige„(7r) structure of the lowest in en- 
ergy ^Eu term of the two-hole Cu04~ center, whose sep- 
aration from the ground state ZR singlet is only about 
1.5 2.0 eV as observed in photoemission experiments^ 
for Sr2Cu02Cl2. Thus, one should expect rather moder- 
ate intensities for the appropriate excitonic mode. On the 
other hand, one might naturally expect a rather strong 
high-energy transition to the predominantly fe°j^e„((T) 
state. Overall, the k dependence of the transition matrix 
elements for the 6jge„-channel is rather typical for purely 
oxygen contributions: the dipole-allowed P modes gradu- 
ally loose their intensity up to zero going from the F-point 
to the BZ boundary both in (tt, tt) and (tt, 0) directions, 
contrary to the dipole-forbidden S modes, whose inten- 
sity rises by approaching the BZ boundary. 

It should be noted that the relative energy position 
of the four CT excitations in the 6ige„-channel is deter- 
mined by the euin) — eu{cr) and the b\g — big separations, 
which are of the order 5.0 and 6.0 eV, respectively, as 
predicted by the above quantum-chemical cluster calcu- 
lations. 



V. DYNAMICS AND DISPERSION OF SMALL 
EXCITONS 

A. Rotational and translational motion of 
two-center excitons 

In contrast to one-center Frenkel excitons, the motion 
of its two-center counterparts is more complicated. We 
addressed above the internal electron-hole motion result- 
ing in 5 — P splitting. Now, we consider the rotational 
and translational motion of two-center excitons. The 
elongated structure of the two-center exciton results in a 
set of both rotational and translational modes, depend- 
ing on the dimension of the lattice of CUO4 centers. Such 
modes for the 2D lattice were addressed in the frames of 
the simple ZN-model»i2i2& these are rotations of the hole 
(electron) around the electron (hole) by 90° and 180°, 
and axial translations (see Fig. 0J). Strictly speaking, 
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FIG. 4: Different types of next-nearest (nn) neighbor or sec- 
ond nearest neighbor transfer, as well as next-next nearest 
(nnn) neighbor or 3rd nearest neighbor transfer of two-center 
excitons. 

Zhang and Ng in their simple model considered only the 
90° and 180° rotation around the electron center with 
the respective transfer integrals ti and ^4, and the ax- 
ial translation accompanied by a 90° rotation around the 
electron center with transfer integral t2- 

It is interesting to note that the 90°- and 180°-rotation 
of the hole around the electron corresponds to the 2nd 
nearest neighbor and 3rd nearest neighbor hopping of 
the hole Cu04~ center which is coupled to a correspond- 
ing motion of the electron CuOj" center. This prob- 
lem is closely related to the well-known problem of the 
ZR-singlet being the ground state of the Cu04~ center, 
moving in the lattice formed by the Cu04~ centersi^ 
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1. Exciton transfer integrals 

Let us consider a three center system ABC and let 
us compare the collinear 180° geometry with the rect- 
angular 90° one. The exciton transfer integrals between 
two-center excitons centered at AB or BC, respectively, 
can be expressed through matrix elements of the effec- 
tive Hamiltonian -ffe// that incorporates potential and 
kinetic energy contributions: 

for diagonal SS and PP transfer, and 

±l(y(.3)_yp))+ (33) 

for off-diagonal SP transfer. Here, are the elec- 

tron (hole) transfer integrals. One has to distinguish the 
collinear geometry with an exciton transfer to 3rd neigh- 
bors Tj^^jl separated by ~ 4A, from the rectangular ge- 

(2) 

ometry with a transfer to 2nd nearest neighbors T^^ 
separated by « we neglected higher 

order terms which are less important. 

According to Eqs. we introduce a set of trans- 

fer parameters to describe the exciton dynamics: 

T5P-i(ri3^-Tf), 
for the collinear exciton motion, and 

rp^V ^ rpXy ^ 1 lrp(2) . rp{2)^ 

^sp ~ '» ' ' 

corresponding to a 90° rotation of the exciton. 

All these parameters have a rather clear physical 
sense. Moreover, the electron (hole) transfer integrals for 

collinear exciton transfer T; are believed to be smaller 

(2) 

than Tg ^ integrals for rectangular transfer. In other 
words, the two-center excitons prefer to move " crab-like" , 
rather than in the usual collinear mode. This implies a 
large difference for the excitonic dispersion in [10] and 



[11] directions. The electronic wave function in the ex- 
citon (contrary to the hole one) has a dominant Cu 3d 

nature that implies a smaller magnitude of the Ti^'"^^ pa- 
rameters compared to the T^'^'^ ones. 

As it was shown by Zhang and Ng,^^'^* the singlet two- 
center S, P exciton can move through the antiferromag- 
netic lattice rather freely in contrast with the single-hole 
motion. Nevertheless, their values for the second order 
integrals ^3 « 0.85 eV and ^4 « 0.65 eV as derived from 
an analysis of the experimental EELS data in terms of the 
simple ZN-model, might be considerably overestimated. 
As will be seen in the more realistic analysis of EELS 
data presented below, there are several channels which 
contribute to the exciton dispersion, such that the nu- 
merical values of each of the exciton transfer parameters 
might be substantially smaller than the effective transfer 
integrals ii^4 of the simple ZN-model. 

Up to now, we addressed the two-center exciton 
dynamics in corner-shared systems like Sr2Cu03 or 
Sr2Cu02Cl2. However, in ID systems such as SrCu02 
there is an additional type of nearest-neighbor transfer 
of two-center excitons localized on the CU2O7 cluster (as 
presented in Fig.^. 



2. Tetragonal symmetry and two-center excitons in the 2D 
Cu02 layer 

The tetragonal symmetry of the 2D Cu02 layer leads 
to certain modifications of the results obtained above 
for excitons in the simple CU2O7 cluster. There, one 
had to distinguish S- and P-like two-center excitons cen- 
tered at the common oxygen site, in analogy to the sym- 
metry properties of O 2s and O 2pa orbitals, respec- 
tively. In contrast to that, two-center excitons in the 
2D Cu02 layer could be readily classified with regard to 
irreducible representations of the tetragonal D^h group. 
Thus, we can form a set of "large orbital" excitons with 
^ig, Big and Eu symmetry (or equivalently with 5, D 
and Px^y symmetry, respectively) rotating around the Cu 
site (with CusOig geometry). Each two-center exciton 
classified by the two-hole state F of the CuO|~ center, 
give rise to one set of 4 excitons. 

The symmetry of the excited CT states interacting 
with the ground state has a crucial infiuence on the tran- 
sition matrix elements. So, given the big symmetry of 
the ground hole state in the isolated CUO4 plaquette, it 
is the "large orbital" charge transfer state with the same 
Big symmetry generated by Sx,y two-center transitions 
that can be mixed with the ground state of the large 
CU5O16 cluster. In terms of the 4 small Sx.y,Px.,y exci- 
tons in the 2D elementary cell this corresponds to the 
mode S- = -^{Sx — Sy). Naturally, the large Eu exci- 
tons generated by small Px^y excitons are dipole-allowed, 
while both Aig and Big excitons are dipole-forbidden. 
In terms of the elementary cell modes the dipole-allowed 
ones are the Px,y, or P± = -j^{Px ± Py) modes. 
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The mixing of the S- mode to the ground state in 
the 2D system defines the intensity of the two-center 
exciton in the Cu02 layer, similarly to the simple Sx 
mode in ID systems. The transition matrix elements 
for dipole-allowed excitons in 2D systems is expressed 
through dipole matrix elements {S^\d\Px^y) instead of 
{Sx,y\d\Px^y) in ID systems; in other words they would 
be a factor \/2 smaller. Thus, the expected intensity of 
the dipole-allowed two-center CT transitions in 2D sys- 
tem is factor 2 smaller than in ID systems. 

B. The dynamics of two-center excitons in ID and 
2D cuprates 

In general, we will consider the two-center F- 
excitons, whose dynamics in frame of the Heitler-London 
approximation^ could be described by an effective one- 
particle excitonic Hamiltonian with standard form 

Heccc= J2 Sri(Ri)T^rir.(Ri-R2)^r.(R2) (34) 
rir2RiR2 

in a site representation with Bp_^ (Ri ) / (R2) being the 
excitonic creation/annihilation operators, or 

ffe.c- 4,(k)Tr,:r.(k)Br.(k) (35) 

in k representation. Here the Fi 2 indices label different 
S- or P-excitons. One should note that instead of the 5*- 
P manifold we could consider also the dynamics of large 
orbital excitons with a certain point group symmetry at 
the F-point. 

1. The dynamics of isolated SP excitonic doublet in ID 
cuprates 

We consider the ID cuprate as a linear chain of corner 
shared CUO4 centers directed along the x-axis. The T(k) 



matrix for an isolated doublet of Sx = S and Px = P 
excitons in the ID cuprate can be written in a rather 
simple form^ 

rp(,^^fEs + 2Tscoskx -2iTspsiiikx \ fof:\ 
^ ' \ 2iTs p sin kx Ep + 2Tp cos kx ) ' ^ ' 

The dispersion of the two-center exciton in ID system is 
governed by the nnn exciton transfer. As it was empha- 
sized above, both {Ts — Tp) and Tsp are believed to be 
small that results in the near-degeneracy of the S and P 
modes with rather trivial dispersion. 

The intensity of the dipole-allowed P exciton is ex- 
pected to decrease sharply with the momentum, while 
that of the dipole-forbidden S'-exciton exhibits a more 
complicated behavior with zeros at the F-point and the 
BZ boundary, and the maximum midway between the 
center and the boundary of the zone at fc = 2/3kmax- 
Overall, one should note that all isolated SP excitonic 
doublets can manifest a rather small dispersion, how- 
ever, the k-dependence of the EELS intensities could be 
extremely large, and they differ significantly for different 
transitions. Naturally, in practice one should account for 
inter-exciton coupling. It is interesting to note that the 
S- and P-excitons are coupled by dipole transition with 
large matrix element, which along with the SP degen- 
eracy makes the SP excitonic doublet in ID cuprates a 
very promising system for nonlinear optical devices^^ 



2. Dynamics of isolated SP excitonic quartet in 2D 
cuprates 

In the 2D case of an ideal CUO2 layer we deal with two 
types of X- {Sx,Px) and y- {Sy,Py) oriented S, P excitons 
in every elementary cell. The T(k) matrix for an isolated 
quartet of Sx,y and Px,y excitons in 2D cuprates can be 
written by a slight modification of the simple ZN-form2& 



T(k) = 



rpxy 



Es + 2Ts COS K 
2iTsp sin kx 
(1 -I- a*{kx, fc, 



E, 



)) 



2iTsp sin kx 
3 -I- 2Tp cos kx 



Tg'"{l + a{kx,ky)) 



nUl-b*{kx,ky)) 



rpxy 
^SP 



spy 

\r||,(i + b*{kx,ky)) rp^ii - a*{kx,ky)) 



(1 - b{kx, ky)) 
Eg + 2Ts cos ky 
2iTsp sin ky 



b{kx, ky)) 
fc,,)) 



Tp (1 — Cl{kx, iiy 

-'2iTsp sin ky 
Ep + 2Tp cos k 



(37) 



y 



where a{kx,ky) — e*'^=" -I- e"*'"'", b{kx,ky) — e*'^^ — e^**^". describe its coupling. For the [11] direction the exciton 
Here, the two 2x2 blocks on the diagonal are related to transfer matrix breaks up into two similar 2x2 blocks 
Sx,Px, Sy,Py excitons, respectively; off-diagonal blocks 
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Eg + 2Ts cos K ± + 2 cos k^) 

T{-2iTspsmk^ - T||,(l - 2i sink^)) 



r 



=F(2zTspsin k^ 



T|^(l 



2i sin kx)) 



Ep + 2Tp cos kx ± rp^(l - 2 cos k^ 



(38) 



with basis vectors 

l^±) = -^iiSx) ± \Sy)y, |p±) = -i=(|p.) ± IP,)), 

(39) 

respectively. In other words, we have two hybrid S+P+, 
and two hybrid S-P- modes. Generally speaking, all 
these modes have different energies at k = (the F- 
point). Both hybrid S-P- excitons involve the dipole- 
allowed (for k || [11]) P_ mode and can manifest itself in 
optical spectra. 

In a nearest-neighbor approximation for the exciton 
transfer, when T5 = Tp = Tgp = 0, Eg ^ Ep, r|| = 
Tpp = T; Tgp = Ti we obtain four modes with energies: 

E^ = E±T +[4. cos^ kx{T^ - T^) + ^T^] * ; 

E^ ^ E±T -[4, cos^ kx (T^ - Pf ) + SPf ] * 
and eigenvectors: 



\{SP) 



s"^ cosalS-i) 



sina|P±); 



with equal weight of S and P modes in eigenvectors. 

In frame of the nn-model for the exciton transfer along 
the [10] direction we address the only example with zero 
SP-mixing when the exciton dynamics is governed by 
the only nonzero parameter T. Due to the "crab-like" 
motion we come to four twisting modes with Sx — Sy and 
Px — Py mixing and energy 

P5± = Es±T[{2+coskxf+smHx)]^- Ep^ = Ep±\T\. 

Here, the dipole-allowed P± modes appear to be disper- 
sionless. Interestingly to note, that instead of one ex- 
pected dipole-allowed mode {Px, or Py) along this direc- 
tion, we obtain two dipole-allowed modes due to Px — Py 
mixing. Naturally, above we addressed rather simple 
"toy" models, however, they allow to clearly demonstrate 
different probable scenarios for the exciton dynamics and 
its manifestation in optical and EELS spectra. 



3. The interaction of two-center excitons in 2D cuprates: 
The interference effects 



\{SP) 



s'"^ sina|5±) 



cosa|P±), 



where 



ri[l + 4sin2 kx]^- 
tan2a= ; ; tant 

21 COSKr 



= 2 sin kx ■ 



It should be noted that the complex matrix elements in 
T± given kx = ky ^ evidence the appearance of " twist- 
ing'^ excitonic modes. If we can neglect the nnn electron 
transfer in comparison with the hole transfer, then we ob- 
tain T = Ti and unexpectedly four dispersionless modes 
with energies: 

Ef ^E±T + V5T; Ef ^E±T - VhT, 

though the eigenvectors continue to depend on kx value. 

If the electron and hole transfer integrals have the same 
magnitude (Te — Tj^ ), then Ti = and we obtain pure 
\S±, P±) modes with a rather conventional dispersion: 



P 



E±T{l + 2coskx); E$ ^ E ±T{1 - 2coskx 



(2) (2) 

Finally, we address the situation with Te = — , 
when T — 0, and we have two doubly degenerate twisting 
modes 

= E-hTi[l+4:Sm^ kx]^; E^ = P-Ti[l+4sin^ /c^^]^ 



Any two excitons of the same symmetry interact with 
each other. The structure of the interaction matrix for 
two Sx^y,Px.y quartets is similar to (|37|l . The inter- 
exciton coupling effects its dispersion, and the intensity. 
The latter is of particular importance for relatively weak 
transitions in neighborhood of the strong ones. Depend- 
ing on the phase of the coupling matrix element, the in- 
tensity of the definite CT transition can either increase, 
or decrease due to its interaction with neighbors up to 
overall suppression. We would like to argue that the par- 
ticular structure of the interaction matrix elements for 
the dipole-allowed Px^y excitons makes the appearance of 
the intensity compensation points a general rule rather 
than an exception. Indeed, the Px^y — Px y coupling is 
governed by matrix elements (Pa,,, |T(k)|P^ y) oc coskx^y 
that change sign midway between F-point and the BZ 
boundary. The Px^y — Py^x coupling along the [11] direc- 
tion is governed by matrix elements (Pa;^j,|P(k)|Py oc 
(1 — 2cosfc2;,j,) that as well change sign at kx,y — \kmax- 

To illustrate inter-exciton coupling effects we address 
two dipole-allowed two-center excitons, say h\g and 6ige„ 

type along the [11] direction with the only nonzero P—P 
coupling. For simplicity, we assume a conventional bare 
energy dispersion for both P modes: 



E{k) = Pp(0) -f Tpcos 



V2 
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with Ep{0) = 3.1, Tp' = 0.4 eV and Sp'(O) = 4.2, Tp = 
0.2 eV, respectively, and conventional EELS intensity dis- 
persion: 



Ipik) = /p(0)- 



2sin^ 



2V2 



fc2 



(1 



/p,(fc)=/p,(0)- 



sm 



2V2 



2fc2 

typical for the big and 6ige„-channels, respectively, with 
intensity ratio like: /p(0) : lp'{0) = 1 : 10. In other 
words, we have neighboring weak low-energy and strong 
high-energy excitons. The interaction matrix element is 
assumed to be 

(p|r(fc)|F') = rpp,(i-cos-^) 

to provide the non-interacting excitons at the F-point. It 
is easy to see that the noticeable interaction effects mani- 
fest itself only in the range {0.5^1.0)kmax- The sign and 
magnitude of Tppi were chosen to provide the intensity 
compensation point near k = 0.7kmax' l^pp'l = 0.3 eV. 
The Fig. |31 presents the model dispersion curves for the 
energy (upper panel) together with those for the EELS 
intensity (lower panel) for two P excitons. We see that 
the interaction results in a rather weak renormalization 
of the bare energy dispersion for both excitons. However, 
the relatively weak interaction leads to a cruciable renor- 
malization of the intensity dispersion for weak low-energy 
exciton with compensation point near q = O.lqmax and 
blazing up near BZ boundary. The appropriate renormal- 
ization for the intensity of the strong high-energy exciton 
is relatively small. Interestingly to notice that the com- 
pensation effect stems from the destructive interference 
of the two excitons and takes place in a rather narrow 
momentum range. In experimental spectra this compen- 
sation can look like an effect of the "lost" band. 

So, we see that the isolated SP excitonic quartets can 
manifest a rather moderate dispersion, however, the k- 
dependence of the EELS intensities is always very strong 
and differs significantly for different transitions. Over- 
all, the model curves in Fig. [3 clearly demonstrate an 
extremely important role played by the EELS intensity 
dispersion. This "matrix element effect", or interference 
effect is likely to complicate the theoretical treatment of 
the experimental EELS and RIXS spectra. 



4- Two-center excitons in 2D systems as compared with ID 
systems 

Concluding this subsection we would like to make a 
short comparative analysis of the main properties of two- 
center excitons in 2D and ID systems with corner-shared 
CUO4 plaquettes: 

1. In ID system we deal only with nnn exciton transfer, 
while in 2D system both nnn and nn transfer exists. 
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FIG. 5: Model dispersion curves for the energy (upper panel) 
and the EELS intensity (lower panel) for two neighboring P 
excitons. For illustration we present the EELS intensity dis- 
persion for isolated S exciton. The momentum dependence 
of the EELS intensity for ftfg-P-exciton demonstrates the de- 
structive interference effect with a full "bleaching" near cer- 
tain critical q-value. See text for details. 



2. At the F-point we have in common 2 dipole- allowed 
transitions associated with one SP excitonic manifold for 
2D system contrary to single transition in ID system. 

3. Integral intensity of dipole- allowed transitions as- 
sociated with one SP excitonic manifold for 2D system 
two times smaller than in ID system. 

4. The S — P splitting in ID system is expected to be 
substantially smaller than in 2D systems. In addition, 
the S — P dipole matrix elements in ID system are ex- 
pected to be as a minimum two times bigger than in 2D 
systems. Both these quantities essentially determine the 
magnitude of the third order nonlinear susceptibility x'^' 
and a potential of the system as an optoelectronic mate- 
rial with high performance Hence, ID systems should 
be more effective nonlinear systems as compared with 2D 
counterparts. Experimental data for ID Sr2Cu03 and 
2D Sr2Cu02Cl2, and exact diagonalization technique on 
small clusters confirm this conclusioniiSiii 
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VI. CT EXCITONS IN OD, ID, AND 2D 
INSULATING CUPRATES PROBED BY EELS 
AND OPTICAL SPECTROSCOPY 

By addressing the nature of the electron-hole excita- 
tions we should emphasize that a comparative analysis of 
the optical and EELS spectra in cuprates with different 
dimensionality of the CUO4 network is extremely inter- 
esting and informative in many respects. In particular, 
the considerable general interest in comparing the differ- 
ent physical properties of the OD, ID and 2D members of 
the growing cuprate family is warmed up by challenging 
problems such as the low-dimensional aspect of electronic 
structure and high- Tc superconductivity. 

The OD cuprates imply a system with well isolated or 
weakly coupled CUO4 plaquettes. Such a situation seems 
to be realized in CUB2O4 which belongs to the class of 
noncentrosymmetric magnetically ordered materials for 
which the unusual coexistence of weak Dzyaloshinskii- 
Moriya type ferromagnetism and inhomogeneous (incom- 
mensurate) magnetic ordering attracts a lot of atten- 
tion. CUB2O4 crystallizes in the tetragonal space group 
IA2dm^ The Cu^+ ions at 4b sites are surrounded by four 
oxygen atoms in planar quadratic coordination so that 
the local symmetry is 4. The Cu^+ ions at 8d sites oc- 
cupy distorted octahedral positions with an exceptionally 
large separation of 3.069 A from the two apical Cu^"*" ions. 
Its local symmetry is 2. 

Many insulating cuprates represent text-book low- 
dimensional model systems. The ID corner-shared or- 
dering of the CUO4 plaquettes is realized in insulat- 
ing Sr2Cu03 to be the best known model system for a 
spin-1/2 antiferromagnetic Heisenberg chain and a very 
promising material for nonlinear optical devices^ 

One of the best representatives for an insulating 
layered 2D copper oxide might be the oxychloride 
Sr2Cu02Cl2, which is iso-structural to the famous 214 
system La2Cu04 but contrary to the latter has a single 
type both of copper and oxygen sites. The compound 
Sr2Cu02Cl2 is one of the most popular model system 
for the insulating phase of the high- Tc cuprates and 2D 
spin-1/2 Heisenberg antiferromagnets, and intensively 
studied both experimentally and theoretically. In this 
tetragonal antiferromagnet (Tn « 250 K) with nearly 
ideal Cu02 planes there are chlorine atoms instead of 
apex oxygens with considerably larger Cu-Clapex separa- 
tion (2.86A) than that of Cu-Oapex (2.42A) in La2Cu04. 
Hence, in Sr2Cu02Cl2 one has the opportunity to ex- 
amine the Cu02 plane states, both copper and oxygen, 
without the "parasitic" contribution of apex oxygens. At 
present, there are a rather large number of experimen- 
tal data for Sr2Cu02Cl2 obtained with the help of op- 
tical spectroscopy,ii X-ray photoemission (XPS),^* ul- 
traviolet photoemission (UPS),^^ X-ray absorption spec- 
troscopy, (XAS)^SiS and angle-resolved photoemission 
spectroscopyi^^i^Si^ML^ 

The OD cuprates with isolated CUO4 plaquettes like 
Bi2Cu04, the corner-shared and edge-shared ID chain 



cuprates like Sr2Cu03 and Li2Cu02, respectively, as well 
as the 2D systems like Sr2Cu02Cl2 should reveal similar 
signatures of the one-center excitations. However, what 
concerns the two-center excitations, there appears a prin- 
cipal difference. Naturally, in tetragonal 2D system of 
corner-shared CUO4 centers, the two-center excitations 
manifest itself equally both in a and b in-plane polar- 
izations, while in ID system of the corner-shared CUO4 
centers they are visible in a similar manner only for " lon- 
gitudinal" polarization when the electric field is parallel 
to the chain direction. For the electric field parallel to the 
CUO4 plane but perpendicular to the chain direction one 
might observe only 7r-like one-center excitations whose 
intensity is determined by the overlap contribution only. 
Thus, by comparing the optical response for an idealized 
ID chain system with two polarizations we could esti- 
mate the contribution of the (T-like two-center excitations 
by simple subtraction of the two spectra. Naturally, it 
would be an oversimplification to say that one can deduce 
the optical response for an idealized 2D plane by simple 
combination of the ID spectra with different polariza- 
tion. A correct extension of the chain data to predict 
the 2D response should also take into account the addi- 
tional 90°-coupling of the two-center excitations already 
at the F-point and the more complex dispersion relations. 
Nevertheless, a simple naive comparative analysis of the 
chain and plane optical responses could provide impor- 
tant semi-quantitative information on the relative con- 
tribution of one- and two-center excitons, their spectral 
positions, and the relative magnitudes of the overlap and 
covalent contributions to the dipole matrix elements. 

One should note the provisional nature of a certain di- 
mensionality for real cuprates. Even in the so-called OD 
systems like CUB2O4, the CUO4 plaquettes are not really 
isolated. It might be surprising that among the best can- 
didates for OD cuprates concerning the optical response 
are such systems like Li2Cu02 or CuGeOs. Usually they 
are considered as ID systems with edge-sharing CUO4 
plaquettes according to their crystalline structure. For 
instance, at room temperature, CuGeOa belongs to the 
orthorhombic group Dg/j (Pbmm). The unit cell con- 
tains two edge-sharing strongly deformed CuOe octahe- 
dra, with two types of Cu-0 bonds of 2.77 and 1.94 A, 
forming one-dimensional antiferromagnetic Cu-0 chains 
along the c axis. However, since the Cu-O-Cu bond an- 
gle along the chain is almost 90°, the transfer from one 
plaquette to its nearest neighbors is strongly suppressed 
for the in-plane states. Thus, the in-plane charge exci- 
tations are localized within one plaquette and with re- 
spect to the electronic properties the CUO4 centers in 
these compounds can be addressed to be approximately 
isolatedpS This standpoint is fairly well confirmed in op- 
tical and EELS data for Li2Cu02i^2*^ In what concerns 
the perpendicular out-of-plane O 2pz orbitals these sys- 
tems should be considered as typical ID chain systems 
with O 2pjr bonding. 

In addition, one should note the problem related to the 
different symmetry of the crystalline field for the CUO4 
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plaquettes in OD and 2D systems as compared with ID 
systems, and between edge-shared and corner-shared ID 
systems. For instance, two inequivalent oxygen sites are 
clearly seen in the polarization-dependent 1s x-ray ab- 
sorption spectra for the ID cuprate Sr2Cu03.®'^ How- 
ever, the energy separation is found to be as small as 
0.5 eV, whic is sufficiently smaller than values of the or- 
der 1.5 — 2.0 eV, expected from predictions of the simple 
point-charge model, or LDA calculations This experi- 
mental finding evidences puzzlingly small non-tetragonal 
effects near the ground state for the CUO4 centers in 
the ID cuprate Sr2Cu03 and allows to make a compara- 
tive analysis of the localized excitations with those of the 
tetragonal 2D cuprate Sr2Cu02Cl2. It seems that the in- 
ternal Cu-0 covalent bond for a Cu04 plaquette provides 
the main contribution to its electronic structure. 



A. The main predictions of the model theory 

Before addressing the experimental EELS data we 
would like to summarize shortly the main results of the 
model theory of one- and two-center CT excitons in 2D 
insulating cuprates. 

1. One-center excitons 

i) We predict three types of dipole-allowed one-center 
excitons: one exciton 6211 with out-of-plane polarization, 
and two excitons eu(7r) and e,i((T) with in-plane polar- 
ization and predominantly vr and cr O 2p orbital weight, 
respectively. 

ii) The relative value of transition matrix element for 
the low-energy (tt) exciton with an estimated energy « 
2 eV and for the high-energy (cr) exciton is determined 
by the level of the vr — cr mixing. 

iii) The low-energy (near 2 eV) part of the absorption 
spectrum for a single CUO4 plaquette is formed by an in- 
terplay of forbidden d—d transitions {h\g h\g, a\g,e^g), 
forbidden b\g — > a2g and allowed h\g — > e„(7r) CT tran- 
sitions, respectively, which are all close in energy. It 
is worth noting that for cuprates, as for many other 
strongly covalent oxides, the conventional division of op- 
tical transitions into crystal-field d—d and CT p—d tran- 
sitions becomes questionable, since the h\g h\g, a\g, 
transitions are accompanied by a strong p—d and p — p 
charge transfer. The low-energy absorption band is ex- 
pected to be a result of strong electron-vibrational cou- 
pling which lifts the selection rules for certain transitions. 

iv) The final eu(7r) and e„(cr) states are unstable with 
respect to the formation of Jahn- Teller (or pseudo-Jahn- 
Teller) centers with localization of the appropriate CT 
excitation. This effect is especially important for the 
b\g — > en(7r) excitation, because the predominant tt na- 
ture of the hole state results in small transfer integrals 
and a relatively large effective mass for the exciton. The 
localization is accompanied by the vibronic reduction of 



transition matrix elements. The high-energy eu{(j) exci- 
ton is expected to manifest significant intensity both in 
optical and EELS spectra. 



2. Two-center excitons 

i) All the two-center excitons generated by the CT 
governed by the strongest cr bonds can be divided into 
three channels foigOig, 6ige„, respectively, in accor- 
dance with the symmetry of the final two-hole state. The 
main effects in optics and EELS are associated with the 
big- and &ige„-channels. In the frame of each channel 
we deal with different final two-hole states whose energy 
and wave functions are mainly determined by correla- 
tions, or configuration interaction effects. The different 
theoretical model calculations predict the lowest energy 
of two-center CT transition for the fe^^-channel to be of 
the order 2 3 eV. 

ii) All the two-center excitons could be divided to even 
S- and odd P-excitons. For 2D cuprates these exci- 
tons generate two sets of four Aig, Big, plane modes 
{S,D,P^^y modes) 

iii) The transition matrix elements in optics and EELS 
have rather complicated k dependence, substantially dif- 
fering for the oxygen and copper contributions. The sim- 
ple dipole approximation fails to explain correctly the 
k-dependent effects in EELS. 

iv) In frame of the big-channel we predict two most 
important excitons with an energy separation « 7.0 eV. 
The low-energy exciton with ZR-singlet as a final hole 
state has rather unusual k dependence of intensity along 
the (tt, tt) direction with a sharp decrease, or even com- 
pensation point, by going to the BZ boundary. The 
high-energy exciton with final ZR-singlet- like ^Aig state 
of predominantly O 2p nature is likely to manifest the 
largest intensity for the &fg-channel. 

v) In frame of the &ige„-channel we predict four most 
important excitons, generated by CT to eu(7r) and e„(cr) 
orbitals with smaller and larger weight of O 2pcr orbitals, 
respectively. These weights determine the relative in- 
tensities of the appropriate excitons. The photoemission 
data^^i^ predict the energy separation between the low- 
lying b\g —> e„(7r) and b\g — > big transitions to be of the 
order 1.5 2.0 eV. The 6ige„ excitons manifest a rather 
simple dispersion of intensity along both main directions 
in BZ zone with gradual interchange between dipole al- 
lowed and forbidden modes. 



3. CT excitons and magnetic subsystem 

All parent ID and 2D copper oxides with corner-shared 
Cu04 plaquettes are antiferromagnets with spin- 1/2, lo- 
calized on Cu atoms of the chains or planes, with ex- 
change integral / w 0.1 eV and a Neel temperature 
Tjv ^ 300 K in 2D. So, we have to address the excitons 
in antiferromagnets. Historically, this problem has been 
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intensively studied in connection with forbidden d—d ex- 
citonic transitions in weakly covalent antiferromagnetic 
insulators like MnF2 (see e.g. the paper by Tanabe and 
Aoyagi in Ref.'s^. At first glance, the main allowed elec- 
tric dipole transitions preserve the spin state and cannot 
result in conventional or orientational spin fluctuations. 
On the other hand, the charge transfer transitions are ac- 
companied by a strong spatial redistribution of spin den- 
sity, and, in a sense, result in a spin density fluctuations. 
Indeed, the two-center CT exciton is formally associated 
with spin singlet-singlet transition from a spin-singlet ini- 
tial state with a nonzero spin density on both centers 
to a spin-singlet final state with a zero spin density on 
both centers. In a sense, this transition is accompanied 
by a local annihilation of the Neel-like or spin-dimerized 
state, and the creation of a non-magnetic vacancy. Such 
a strong spin density fluctuation could result in effec- 
tive two-magnon (2M) processes. In other words, one 
might expect strong 2M Raman scattering generated by 
a two-center exciton. Indeed, the experimental Raman 
measurements for different insulating cupratesifl^ show 
that the maximal strength of 2M scattering occurs for 
excitation energies substantially higher than the optical 
band edge with a peak which we assign to a two-center 
exciton. 

For small one-center excitons we have a rather conven- 
tional s = l/2— >s = l/2 transition with spin density 
fluctuation localized inside the CUO4 plaquette. Such a 
transition is not accompanied by strong 2M Raman pro- 
cesses, what could be used to identify such kind of exci- 
tons. The redistribution of spin density from the copper 
atom to the oxygen ones for the big e„ transition 
switches on the strong ferromagnetic Heisenberg O 2p- 
Cu 3d exchange with the nearest neighbor CUO4 plaque- 
ttes. Interestingly, the different sign of exchange coupling 
for e„ and big holes with the same neighborhood, ferro- 
magnetic for the former, and antiferromagnetic for the 
latter, leads to a number of temperature anomalies near 
the 2D-3D antiferromagnetic phase transition. There is 
at first the blue shift effect for the energy of the small 
one-center exciton big by lowering the temperature 

near and below Tjy. Indeed, at T ^ Tjv, the average 
molecular field for the CUO4 center turns to zero. The 
3D antiferromagnetic ordering is accompanied by a rise 
of exchange molecular fields and respective spin split- 
tings. Due to the different signs of molecular fields for 
e„ and big states this is accompanied by an increase of 
the transition energy with maximal value of blue shift 
A « /3(|iJf,jg| -I- |i?e„|). This quantity could be as large 
as several tenths of eV. Additionally, one has to expect 
a strong (of the same order of magnitude) broadening of 
the excitonic line with increase of the temperature due 
to strong fluctuations of molecular fields. All these ex- 
pectations are experimentally found for the 2.0 eV line in 
Sr2Cu02Cl2^ confirming its one-center excitonic nature. 
A similar situation is observed in La2Cu04)^ although 
the authors have explained the data by assuming a po- 
laronic nature of electrons and holes with a short-range 



interaction in between. 

It should be noted that at present there is no relevant 
theory of spin-excitonic coupling for small CT excitons 
in strongly covalent and correlated systems. The tra- 
ditional framework for understanding the two-magnon 
Raman scattering in antiferromagnets for non-resonant 
excitation energies has been an effective Loudon-Fleury 
Hamiltonian,^^ which implies well localized spins and 
weak spin-wave like spin fluctuations. 



4- CT excitons and phonon subsystem 

Generally speaking, the CT exciton creation is usually 
accompanied by an excitation of lattice modes. Indeed, 
the hole transfer from Cu 3d to O 2p state, or from an 
ionic to a covalent configuration is accompanied by a sig- 
nificant shortening of the equilibrium Cu-0 bond length. 
The exciton-phonon interaction strongly influences the 
line-shape of absorption and results in a phonon Raman 
scattering. The measurement of the Raman intensity as 
a function of excitation light energy is a very informative 
probe of the origin of electronic transitions.^" In particu- 
lar, this method has allowec^ to resolve a fine structure of 
the low-energy (LE) excitonic feature in La2Cu04 with 
a sharp peak at 2.14 eV as narrow as 50 meV and a 
broader structure at 1.9 eV. In addition, the authors 
made the important observation of the selective reso- 
nance enhancement, namely that neither the first-order 
even-parity phonon nor their harmonics show resonant 
enhancement near the LE exciton peak, while the high- 
energy (above 450 cm~^) presumably odd Raman modes 
show large resonant enhancement of as much as a factor 
of 10 — 40 relative to a laser energy of 2.7 eV. Similar re- 
sults have been obtained for insulating YBa2Cu306.if^ 
In our opinion, such an unusual behavior of the LE ex- 
citon could be associated with its Jahn- Teller nature. 
Indeed, the big —^ e„ transition to the orbital doublet 
Cu represents a textbook example of a so-called A — E 
transition.^* The excited, orbitally degenerate e,i state is 
unstable with regard to vibronic coupling with the local 
distortion modes Aig, Big, B2g, and the formation of 
a polaron-like (soliton-like) vibronic center with a com- 
plex two- or four-well adiabatic potential and a rather 
strong renormalization of vibration frequencies. In other 
words, the doublet e„ state tends to a spontaneous lo- 
cal symmetry breaking, including removal of the inver- 
sion center due to interaction with the close in energy 
even states and pseudo- Jahn- Teller effect. Namely the 
latter would result in a strong resonance coupling of the 
big — > e„ exciton with odd E^ lattice modes. Naturally, 
the structure of such a center would strongly depend on 
differences in the bare lattice and elastic parameters and 
differ in 214 and 123 systems. One should emphasize 
that the formation of the heavy polaron-like, or localized 
small exciton would result in a strong enhancement of its 
effective mass. 
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B. Cuprates with well isolated or weakly coupled 
Cu04 plaquettes (OD cuprates) 
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FIG. 6: Low-energy absorption spectra for the OD cuprate 
CuB204(k II z,E _L z)^^. 

In Fig.Elwe reproduce the low-energy absorption spec- 
tra for the OD cuprate CuB204.*^^ In contrast to most ID 
and 2D cuprates where the fow-energy absorption bands 
are broad and featureless, the absorption spectra for the 
true OD cuprate CUB2O4 (Fig.O provide the remarkable 
opportunity to observe the interplay of low-energy al- 
lowed and forbidden transitions. The experimental spec- 
trum shows a number of narrow zero-phonon peaks with 
up to 70 well-resolved phonon sidebands^' imposed on 
rather broad bands. The former may be assigned to for- 
bidden b\g — > b2g, a\g,e^g {d - d) and b\g a2g {p - d) 
transitions, while the latter may be assigned to allowed 
b\g Bu (tt) CT transitions in two types of CUO4 plaque- 
ttes, respectively. The transition energies at Cu4b04 and 
Cu8d04 plaquettes are similar because of the small in- 
fluence of the remote apical ions. The experimental 
EELS spectra for CUB2O4™ exhibit a rather broad dis- 
persionless band peaked near 5 eV which could be associ- 
ated with a phonon-assisted localized electron excitation 
b\g e^. The decrease of intensity with the increase 
of momentum agrees with the cluster (CUO4) model for 
electronic states and the dipole nature of transition. A 
weak feature distinctly observed at small momenta near 
2 3 eV could be assigned to the localized electron ex- 
citation b\g — > e^, i.e. to the low-energy bonding coun- 
terpart of e„ state. The optical and EELS spectra for 
OD cuprates are governed only by the intra-center tran- 
sitions. 

Interestingly, that very similar EELS spectra are ob- 
served in Sr2Cu03 for polarization q || b, that is perpen- 
dicular to the chain direction {out- off- chain spectra;2i 
Fig. [71). Indeed, like for OD systems, there are no 
two-center e/i-excitations governed by strong covalent a- 
bonding in this case, and we deal with the predominant 
contribution of intra-ccnter excitations. So, for the polar- 
ization perpendicular to chain direction, the ID cuprates 
like Sr2Cu03 provide the optical and EELS response typ- 
ical for a OD system. Contrary to CuBi204, both dipole- 
allowed b\g ef;° transitions in Sr2Cu03 manifest itself 
more distinctly with well-defined EELS-peaks at 2.0 and 



5.4 eV, respectively. The appropriate peaks in optical 
conductivity are situated at 1.8 and 4.3 eV, respectively. 
One should note that these EELS data present a straight- 
forward experimental manifestation of the dipole-allowed 
intra-center b\g e^'" excitations without the "para- 
sitic" effect of inter-center transitions. 

The situation in ID systems crucially changes for po- 
larization parallel to the chain direction, or for the in- 
chain spectra. Even in Li2Cu02, where, similarly to real 
OD systems, we deal with a predominant contribution 
of intra-center transitions, the in-chain EELS spectra— 
look substantially different from those for CuBi2 04. First 
of all, it concerns the line-shape of the intensive band in 
the spectral range 4-^ 5 eV with a very sharp peak at 4.7 
eV (F-point) and small negative dispersion, which evi- 
dences its excitonic nature. Indeed, the concept of the 
one-center b\g e° exciton was successfully appliedS^ 
for a quantitative description of these EELS spectra. 

Similar absorption spectra as shown in Fig. are also 
observed in ID cuprates with 90° Cu-O-Cu bonds like 
CuGeOafS? The EELS spectra of Li2Cu02 present an- 
other example for the weak energy dispersion of the one- 
center exciton accompanied by a strong dispersion of its 
intensity. The different behavior of the one-center dipole- 
active b\g —> e„ excitations with in-chain and out-off- 
chain polarization is straightforwardly associated with 
specific translational symmetry and lattice dynamics of 
ID systems. Indeed, only in-chain component of mo- 
mentum is conserved and can describe the electron and 
lattice modes. Additionally, the ID systems are soft with 
respect to out-off-chain distortions that favor the local- 
ization of electron excitations. The most probable mech- 
anism governing the structure of the polaron-like one- 
center exciton in Li2Cu02 is associated with strong cou- 
pling of the purely oxygen O 2p„ hole mode and oxygen 
distortion modes at the BZ boundary. 



C. Polarization dependent EELS spectroscopy of 
ID copper oxide Sr2Cu03 and separation of one- 
and two-center CT excitons. 

In 2D systems we usually deal with spectra being a 
hardly resolved superposition of both types of excitons. 
This renders the separation of transitions with differ- 
ent nature rather difficult. Therefore, there remains still 
some ambiguity concerning the reliable identification of 
two dipole-allowed one-center CT excitons and seven es- 
sential two-center CT excitons that might still question 
its existence as well-defined entities. Unfortunately, this 
concerns also the structure of the low-energy optical re- 
sponse observed in the spectral range 2-;- 3 eV, which is of 
special importance since it is associated with the states 
which are believed to define the unconventional proper- 
ties of the cuprates. One should emphasize that conven- 
tional optical measurements are restricted only to the F- 
point and cannot separate between localized one-center, 
dispersionless excitations and two-center excitons with a 
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noticeable dispersion. In contrast to optics, the angle- 
resolved EELS spectroscopy provides unique opportuni- 
ties to reveal the exciton dispersion and separate one- and 
two-center CT excitons. One-dimensional cuprate com- 
pounds are good candidates for such a study since the 
scattering of electrons with a transferred momentum q 
perpendicular to the chain direction excites only electron- 
hole pairs sitting on one CUO4 plaquette. On the other 
hand, for q parallel to the chain direction, both types 
of excitons (one- and two-center) can be observed. The 
interpretation of the EELS spectra with "longitudinal" 
response were already proposed in terms of the stan- 
dard one-band Hubbard modeli22i Below we show that 
the polarization-dependent angle-resolved EELS study of 
the ID cuprate Sr2Cu03 with corner-shared CUO4 pla- 
quettes provides a unique opportunity to separate both 
the one- and two-center CT excitons, and the two types 
of one-center excitons, with the first unambiguous man- 
ifestation for the relevance of the O 2p^ holes for the 
low-energy excitations in cuprates. 

The EELS spectra for Sr2Cu03 in "longitudinal" re- 
sponse with the transferred momentum oriented along 
the chain direction were measured earlier onS and have 
been interpreted within standard Hubbard models. How- 
ever, such models can describe properly only the "longi- 
tudinal" response with the transferred momentum ori- 
ented along the chain direction. In this paper we argue 
that along with the "longitudinal" response a ID cuprate 
compound like Sr2Cu03 should reveal rather strong low- 
energy " transversal" response for the transferred momen- 
tum oriented in the plane of the CUO4 plaquettes but 
perpendicular to the chain direction. The EELS spec- 
tra for Sr2Cu03 in both polarizations are presented in 
Fig. 13 Their comparison leads to very important con- 
clusions. As is to be expected, we see a strong difference 
between the two sets of spectra. First, this concerns the 
well-defined dispersionless EELS peaks at 2.0 and 5.5 eV 
in transversal polarization (right hand side panel in Fig. 
0). The intensity considerations and the absence of no- 
ticeable energy dispersion allow us to associate them with 
dipole-allowed CT transitions having a particularly local- 
ized nature. Moreover, when comparing both spectra we 
see that despite the strong inequivalence of longitudinal 
and transversal polarizations in Sr2Cu03 the low-energy 
transition peaked near 2 eV is equally present in both 
polarizations. Though for longitudinal polarization this 
is partly hidden for low momentum values by the inten- 
sive band peaked at 2.6 eV and is seen as a shoulder, 
near the BZ boundary it is a well-separated weak band 
due to the big blue-shift of the intense neighbor. This 
fact implies that the excitation is localized on the sin- 
gle Cu04 plaquette being the only common element of 
longitudinal and transversal geometry in this ID cuprate 
with corner-shared CUO4 plaquettes. Hence, by taking 
into account the intensity ratio we can unambiguously 
identify the 2.0 and 5.5 eV peaks in the EELS spectrum 
of Sr2Cu03 with the one-center CT excitons e„(7r) and 
eu{cr), respectively. 
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FIG. 7: EELS spectra in Sr2Cu03 for longitudinal q || a (left 
panel) and transversal q || b (right panel) response with an 
illustration of one- and two-center excitons. 



The sizeable dispersion of the most intense low-lying 
CT exciton peaked in EELS at 2.6 eV agrees with its 



two-center natun 
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that is fairly well confirmed in the 



studies of nonlinear optical effects in Sr2Cu03. In- 
deed, our analysis of the photoinduced absorption ex- 
periments and the measurements of the third-order non- 
linear susceptibility x^"^-* in Sr2Cu03^®'^'* allows to reveal 
the near-degeneracy of both S- and P-types of two-center 
CT excitons with an energy of about 2.0 eV, and to ob- 
tain a reliable estimate for the transition matrix element 
(EH: |(S'|d|P)| w 2e X 4 A (see Fig. 4 in Ref. 18), that 
straightforwardly points to a two-center CT exciton. In 
addition, the nonlinear optical response of Sr2Cu03i^ re- 
veals a weak spectral feature red-shifted to several tenths 
of eV with regard to the position of the main two-center 
exciton. It can be attributed to a weak one-center ex- 
citon that seems yet to be invisible in reflectivity mea- 
surements of Sr2Cu03.^ In this connection, we would like 
to emphasize once more the decisive role of direct EELS 
measurements in the observation and unambiguous as- 
signment of one- and two-center excitons in Sr2Cu03 as 
compared with conventional indirect optical data. 

It is interesting to note that a low-energy transition 
with small dispersion was also seen in the RIXS data of 
the ID cuprate SrCu02 pubhshed in Ref. [s^ In that 
paper it was interpreted as the onset of the continuum of 
the one-band Hubbard model. By comparing it with the 
above discussion of the EELS data of Sr2Cu03 we are 
led to an alternative interpretation of these RIXS data 
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in terms of a superposition of two different electron-hole 
excitations, namely the above discussed one- and two- 
center excitons. 



D. CT excitons in 2D insulating cuprates 

At first sight, all the above findings are restricted to 
ID cuprates. However, already a shorthand inspection 
of the high-resolution EELS spectra for the 2D system 
Sr2Cu02Cl2 for two polarizations22i2i shows nearly the 
same behavior of the lowest in energy electron-hole exci- 
tations. Being encouraged by such an impressive success 
of our model we can now address the detailed analysis of 
CT excitons in 2D insulating cuprates. 

1. EELS spectra of Sr2 CuO^ Ch 

As it was already noted, Sr2Cu02Cl2 is one of the best 
realizations of a 2D antiferromagnetic, insulating model 
compound. Of particular importance for the assignment 
of both one- and two-center excitons are data obtained by 
ARPES which is a powerful tool to examine the energy 
spectrum of the two-hole Cu04~ center. It probes one- 
particle excitations in contrast to optics and EELS, and 
the nonbonding oxygen states are clearly visible. At the 
F-point the selection rules for ARPES are the same as for 
optical and EELS transition, so that ARPES " sees" only 
the purely oxygen photoholes, or the states of 
the two-hole Cu04~ center with 6igei,-like configuration. 
The inspection of the experimental ARPES spectra for 
Sr2Cu02Cl2 obtained in Refs. i43ii44 gives valuable infor- 
mation regarding the two low-lying ^E^ states. It clearly 
reveals two strong bands separated from the ground state 
ZR-singlet by 1.5 2.0 eV and « 5.0 eV, which could be 
naturally related to the low-energy eu(7r) and high-energy 
e„(cr) photohole states, respectively. Different theoretical 
estimations^ii^SiSiZ^iZ^ corroborate the ARPES data and 
point to a rather low-energy position of the O 2-p-p; states 
in ID and 2D cuprates with a corner shared arrangement 
of Cu04 plaquettes with Ao-tt = (cpo- — CpTr) ~ 1 3 eV, 
where t^^ and e^jr are the centers of gravity for different 
O 2pcr and O 2p^ states, respectively. 

Unfortunately, the measurements in Ref. 3 are en- 
ergy restricted and do not provide information regard- 
ing the high-energy 6ige„(7r) and 5"ge„((T) configura- 
tions. However, recent resonant X-ray scattering spec- 
troscopy (RIXS) measurements for different insulating 
cupratcE^^'^"^'^^ with the incident photon energy tuned to 
the Cu K edge (hard X-rays) reveal a rather wide (2-^3 
eV) band peaked at 5 6 eV, which could be unambigu- 
ously identified as the band of antibonding Cu id-O 2p 
(a°g, Bg) statcs. Hence, the b\g — b^g separation 

could be estimated to be « 6.0 eV. 

Thus, making use of the ARPES and RIXS data, as 
well as the theoretical predictions, we can establish an 
overall picture of optical and EELS spectra (at the F- 



point). These spectra are governed by dipole-allowed 
one- and two-center CT excitations and they include: 

i) the lowest in energy (1.5 ~ 2.0 eV) one-center b\g 
et((7r) transition with rather small intensity due to the 
dominantly O 2pt^ nature of the final state and strong 
tendency to self-localization (trapping); 

ii) the high-energy (« 5.0 eV) and rather intensive one- 
center b\g — > e„(cr) transition; 

iii) in frame of the ^^^-channel we predict three two- 
center CT transitions: 1) the rather intensive lowest in 
energy (2.0-;- 3.0 eV) transition with ZR-singlet as a final 
hole state; 2) the most intensive high-energy transition 
with ZR-singlet-like final hole state blue-shifted to « 7.0 
eV; and 3) probably the less intensive transition blue- 
shifted to « 5.0 eV with regard to the first one; 

iv) in frame of the &ige„-channel we predict two dou- 
blets of two-center CT transitions shifted « 6.0 eV with 
regard to each other. The high-energy doublet is rel- 
atively more intensive. The partial transitions in both 
doublets are shifted 3.0 ^ 3.5 eV with regard to each 
other. Among them, the lowest transition which is rather 
intense to the bigCu', dpi: state is blue-shifted by 1.5 2.0 
eV with regard to the lowest in energy big transition with 
ZR-singlet as a final hole state, while the more intensive 
high-energy 6ige„; dpa transition is blue-shifted by w 5.0 
eV with regard to the same transition. In this quartet the 
most intensive CT transition with bigeu',pp(J final state 
is expected to have the maximal energy (~ 13 14) eV 
among all the CT transitions governed by the strong a 
bond. 
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FIG. 8: The spectral dependence of the imaginary part of di- 
electric function e{ijj) for Sr2Cu02Cl2. Arrows mark the pre- 
dicted energy position of one-center, two-center bfg-channel, 
and feigBu-channel CT transitions. 

As we shall see, all these predictions are in agreement 
with experimental optical and EELS spectra available for 
Sr2Cu02Cl2^ 
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In Fig. |S1 we present the spectrum of the imaginary 
part of dielectric permittivity derived from the Kramers- 
Kronig transformation of the EELS dataiii Remarkably, 
all the excitations which were analysed before correspond 
to visible features in e2. The weak low-energy feature 
near 2.2 eV in the EELS spectra (2.0 eV in £2 and optical 
conductivity) could be ascribed to the lowest in energy 
one-center b\g — > eu(7r) transition. This feature occu- 
pies the tail of a rather intensive band peaked at 2.7 eV 
in EELS spectra (2.4 eV in €2 and optical conductiv- 
ity). This band is naturally associated with the lowest 
in energy b'lg big two-center CT transition with ZR- 
singlet as a final hole state. Its energy is of particular 
importance for the whole set of two-center transitions. 
Having positioned the lowest one- and two-center exci- 
tons, the higher states are fixed by the chosen parame- 
ter values (see Fig. EJ- Two features in EELS, near 4.2 
eV and 7.1 eV could be related to two-center bigCu, dpn 
and bigCu'idpa transitions, respectively, while the wide 
band near 6.0 eV in EELS could be naturally assigned 
to the one-center b\g — > e„(tT) transition, whose energy 
is slightly higher than that of its two-hole counterpart. 
The overall spectrum of the most effective CT transitions 
generated by the Cu 3d-0 2p a transfer ends with very 
strong features near 9-;- 10 and 12-^ 13 eV in EELS which 
could be certainly assigned to the high-energy two-center 
b\g',PP CT transition with ZR-singlct-like final hole state 
of predominantly pp configuration, and the transitions 
to the 5ige„;pp7r and bigeu',pp<y final states, respectively. 
All this shows that important spectral information is con- 
tained in the range above 8 eV. The peak at about 18 eV 
is likely to be attributed to transitions with O 2s initial 
state. 

The interpretation of EELS spectra for nonzero mo- 
mentum becomes more complicated due to the energy 
and intensity dispersion of the dipole-allowed modes and 
the appearance of numerous new modes which are forbid- 
den at the F-point. Hereafter, we focus only on angle- 
resolved EELS spectra in the spectral range up to 8.0—8.5 
eV for several momentum values in [100] and [110] direc- 
tions, obtained by Wang et alMi and by Fink et alm^^ 
The spectra differ somewhat by the maximal momenta 
values and the considerably better resolution in the latter 
case. For illustration, we reproduce in Fig. |51 the EELS 
spectra^SiSi along the [110] direction. The low-energy 
part of the EELS spectra in the 2D system Sr2Cu02Cl2 
along the [110] direction and the longitudinal spectra in 
the ID system Sr2Cu03 have a very similar structure up 
to some quantitative coincidence. For both compounds 
the main spectral feature is associated with an intensive 
band related to the two-center b\g ZN-exciton with clear 
dispersion peaked at 2.6 eV near the F-point and at 3.0 
eV near the BZ boundary. In both systems this band has 
the low-energy shoulder near 2.0 eV which is distinctly 
seen in the F-point or near the BZ boundary. This two- 
peak structure of the CT band in Sr2Cu02Cl2 is corrobo- 
rated by conventional optical measurementsiS^i^ (see Fig. 
I10|l . The EELS spectra along the [100] direction manifest 
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Energy loss (eV) 

FIG. 9: The spectral dependence of the loss function for 
Sr2Cu02Cl2. The region of possible destructive exciton- 
exciton interference is highlighted by a circle. See text for 
details. 



an expected overall drop in intensity with relatively small 
energy dispersion. In contrast to the [100] direction, the 
main peak for the [110] direction of Sr2Cu02Cl2 exhibits 
clear signatures of strong dispersion, particularly for mo- 
mentum values in the range {Q.b^Q.7)kmax- Namely this 
effect was a starting point for the model theory of the 
CT excitons by Zhang-Ng.-'^^'^'^ However, in our opinion, 
these signatures are not only related to the strong energy 
dispersion of the ZN-exciton, but also to the unconven- 
tional behavior of the intensity dispersion for different ex- 
citons with the pecularities especially in the momentum 
range (0.5 l.Q)kmax shown in Fig.|31 Indeed, as it was 
emphasized above, the intensity of the main peak in the 
low-energy part of the EELS spectrum for Sr2Cu02Cl2, 
assigned to the dipole-allowed P-exciton associated with 
the b\g big two-center CT transition with ZR-singlet 
as a final hole state, sharply decreases with the increase 
of momentum along [110] direction with a probable com- 
pensation point near the BZ boundary. This circum- 
stance allows to clearly observe a sharp rise of the spec- 
tral weight in a rather broad range with the well-defined 
peak in EELS near 3.8 eV. This spectral feature may 
be unambiguously associated with the dipole-forbidden 
Ag counterpart of two-center bigeu',dpa exciton, which 
is allowed and rather intensive at the (tt, tt) point. In- 
terestingly that the result of this specific behavior of the 
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EELS intensity for different excitons might be mistaken 
for manifestation of the energy dispersion of the main 
peak associated with the g — u dipole-allowed two-center 
big exciton, as it was made by Wang et alM^ NaturaUy, 
this mistake leads to a conclusion about very large (~ 1.5 
eV) energy dispersion for this exciton. The real energy 
dispersion for different excitons in our opinion may not 
exceed the reasonable values of the order of 0.5 ~ 1.0 eV. 

Thus, the overall analysis of the experimental EELS 
spectra in the model 2D insulating cuprate Sr2Cu02Cl2 
allows us to certainly assign a set of distinctly observed 
features in a rather wide spectral range 2.0 -;- 14.0 eV 
to the predicted one- and two-center CT excitons, and 
confirm the validity of the theoretical concept based on 
an embedded Cu04 cluster model. 

Interestingly, that the experimental data available al- 
low us to estimate the difference in the eZ-/i-bonding en- 
ergies Eh for the low-energy one- and two-center CT ex- 
citons. Indeed, 

Ebieuin)) - Ebihgbig) = ECE^) - E^A^g) 

~{E{eu{^)) - E{bigbig)) « (2 ± 0.5) eV, 

where the first term in the right hand side represents the 
energy of big — > e„(7r) excitation in the hole Cu04~ cen- 
ter with ZR ground state, while the second one represents 
the difference in the energies of the respective excitons. 
For the estimate we have made use both of ARPES and 
EELS data. 



2. Fundamental absorption hand in parent cuprates: 
Interplay of one- and two-center CT excitons 

A close examination of other materials shows that the 
two-component structure of the CT gap appears to be a 
common place for all parent cuprates (see e.g., Refs. ITll 
0, and0). The dipole-allowed locahzed big — > excita- 
tion within the Cu04 plaquette related essentially to the 
eu(7r) state is distinctly seen in optical and EELS spec- 
tra for different insulating cuprates as a separate weak 
feature or a low-energy shoulder of the more intensive 
band near 2.5 eV assigned to the inter-center CT tran- 
sition associated with the Zhang-Rice singlet-like excita- 
tion b\g]pd. 

Both components are characterized by a different cou- 
pling to the magnetic and phonon subsystems thus pro- 
viding additional ways to separate them. For small 
one-center excitons we have a rather conventional s = 
1/2 — > s = 1/2 transition with a spin density fluctua- 
tion localized inside the CUO4 plaquette. In contrast to 
the two-center spin-singlet excitation, such a transition 
is not accompanied by strong two-magnon Raman pro- 
cesses, what could be used to identify this kind of exciton. 
The redistribution of spin density from copper to oxygen 
for the big transition switches on a strong ferro- 

magnetic Heisenberg exchange between nearest neighbor 



CUO4 plaquettes. The different sign of exchange coupling 
for the and the big holes (ferromagnetic for the for- 
mer, and antiferromagnetic for the latter) leads to a num- 
ber of temperature anomalies near the antiferromagnetic 
phase transition. First, one expects a line broadening 
and a blue shift for the energy of the one-center exciton 
i>ig ~^ e„ by lowering the temperature near and below 
Tjv. All these expectations are experimentally found for 
the 2.0 eV line in Sr2Cu02Cl2^^ confirming its one-center 
CT excitonic nature. 

The behavior of one- and two-center excitons in 
monoxide CuO^ has some specific features due to the 
low (monoclinic) crystal symmetry. The etj(7r) state 
is splitted into two components with markedly differ- 
ing inter-chain e„(7r) — big exchange, and the respective 
CT bands near 1.7 eV show markedly different tempera- 
ture behaviors?^ Some experimental data which demon- 
strate the two-peak structure of the CT gap in insulating 
cuprates are presented in Fig. EH 



VII. CONCLUSION 

We have developed a semi-quantitative cluster ap- 
proach based on the complete Cu id-0 2p set of the em- 
bedded CUO4 cluster orbitals with a reasonable choice 
of single-particle and correlation parameters to consis- 
tently describe the electron-hole excitations in insulating 
cuprates in a rather wide energy range up to 10 15 
eV. We extended the Zhang-Ng model by considering the 
complete set of Cu Sd and O 2p orbitals and by introduc- 
ing one-center and two-center excitons. In this connec- 
tion, it is worth noting that the small charge transfer ex- 
citon of the ZN-model can actually be attributed namely 
to one-center excitons rather than to two-center ones (see 
e.g.. Fig. 3 in Ref. .28^ . Instead of one d — p^ CT transi- 
tion of the ZN-model with an energy w 2.5 eV (F-point) 
we arrive at a set of one- and two-center excitons gen- 
erated hy d — Per charge transfer, and occupying a very 
broad energy range from « 2 up to « 13 eV. Moreover, 
the largest dipole intensities have the excitons with the 
highest energy 9.0 13.0 eV. 

The simple cluster model allows to consistently ac- 
count for correlation effects in the final two-hole states 
for two-center excitons. It is shown that these effects 
are of particular importance both for the energies and 
intensities of such excitons. 

We argue of an important role played by the transition 
matrix element effects both in optical and EELS spec- 
tra. In the latter case we obtain the momentum depen- 
dence of matrix elements for different both intra-center 
and inter-center transitions. One should note that the 
optical and EELS measurements provide a very impor- 
tant instrument to inspect the electronic structure and 
the energy spectrum of insulating cuprates. However, 
decisive conclusions might be taken only after a detailed 
account for "matrix element" or intensity effects. 

The exciton-exciton interaction is shown to be of par- 
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FIG. 10: Some experimental results of EELS and optical 
measurements for different insulating cuprates which clearly 
demonstrate a structure of the optical gap: a) low-energy 
part of EELS spectra in Sr2Cu02Cl2 for k |j [100}^^; b) op- 
tical conductivity in Sr2Cu02Cl2i"; c), d) imaginary part of 
the dielectric response in Sr2Cu02Cl2iii^and YBa2Cu306," 
respectively, together with the integrated intensity of the re- 
lated two-magnon (2M) Raman spectrum as a function of the 
laser excitation energy; e) imaginary part of the dielectric re- 
sponse in La2Cu04^ and the normalized integrated intensity 
of 1205 cm~^ Raman line in La2Cu04 as a function of the 
laser excitation energy"^ f) contribution of two fitted oscilla- 
tors to the imaginary part of the dielectric function for the 
insulating Eu2Cu04^ 



ticular importance for the actual momentum dependence 
of the EELS intensities being a probable origin of such 
pecuharities as an intensity compensation point due to a 
destructive exciton-exciton interference. 

Starting from the predictions of the model theory 
and some ARPES data we were able to give a semi- 
quantitative description of the experimental EELS spec- 
tra for Sr2Cu02Cl2 that is free from many shortcomings 
of the simplified ZN-model. Both, the different experi- 
mental data and the theoretical analysis show that the 
nature of the CT gap in insulating parent cuprates is 
determined by two nearly degenerate excitations, a lo- 
calized intra-center and a inter-center CT exciton with 
considerable dispersion. The former represents the op- 



tical counterpart of the so-called "1 eV peak" revealed 
by ARPES measurements in a number of cuprates and 
is associated with the hole CT transition big e„(7r) 
from the Cu 3d-0 2p hybrid big oc d,j.2-y2 state to the 
purely oxygen O 2p^ state localized on the CUO4 pla- 
quette, while the latter corresponds to the b'lg big CT 
transition between neighboring plaquettes with the ZR 
singlet being the final two-hole state. 

Making use of the EELS measurements for the ID 
cuprate Sr2Cu03 with two different directions of the 
transferred momentum, we demonstrate straightfor- 
wardly the two-peak nature of the CT gap with the co- 
existence of nearly degenerate one- and two-center CT 
excitons. In contrast to the 2D case, the two different 
types of excitons may be well separated by choosing the 
transferred momentum to be perpendicular to the chain 
direction. 

The structure of the optical gap with two rather well- 
defined CT excitons seems to be typical for a wide group 
of parent cuprates that implies a revisit of some generally 
accepted views on the electronic structure both of ID and 
2D cuprates. 

The nonbonding O 2pt^ states with and a2g sym- 
metry form the energetically lowest purely oxygen hole 
states localized on a CUO4 plaquette with energy « 2 
eV. One might speculate that these O 2p^ states could 
be as preferable for the localization of additive hole as 
the big oc dx2-y2 ground state which would result in the 
instability of the Zhang-Rice singlet in doped cuprates. 
The most probable candidate states for a competition in 
energy with the Zhang-Rice (big) (or ^Aig) singlet are 
the singlet ^Eu, or the triplet ^Eu terms of the bigCu 
configuration. Such a ^Aig — ^'^E^ competition for the 
hole-doped CUO4 plaquette can result in an unconven- 
tional behavior of doped cuprates. Perhaps, this could 
substantially infiuence the nature of the superconduct- 
ing state in these systems. Such a conclusion is sup- 
ported both by local-density-functional calculations^^i^^ 
ab initio unrestricted Hartree-Fock self-consistent field 
MO method (UHF-SCF) for copper-oxygen clusterSfi^iS 
and a large variety of experimental data. To the best 
of our knowledge, one of the first quantitative conclu- 
sions on a competitive role of the hybrid copper-oxygen 
big{dx2_y2) orbital and purely oxygen O 2pTr orbitals in 
the formation of valence states near the Fermi level in the 
Cu02 planes has been made by A.K. McMahan et alJ^ 
and J. Tanaka et alJ^ Namely these orbitals, as they 
state, define the low-energy physics of copper oxides. 

Insulating 2D cuprates like Sr2Cu02Cl2 are not molec- 
ular crystals, and making use of the small exciton ap- 
proach inevitably leads to problems with strong excitonic 
overlap and mixing, and manifestation of band-like ef- 
fects owing to electron-hole pair decoupling. Neverthe- 
less, we suppose that the strong electron-lattice polar- 
ization effects^ for CT states may provide an effective 
localization of CT excitons. So, the small exciton model 
turns out to be a good and instructive approximation to 
describe optical and EELS spectra in different insulating 
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cuprates with CUO4 plaquettes. In this connection, we 
would Hke to emphasize the specific role played by the 
ground state of the cuprate in the formation of optical 
and EELS response. For a strongly correlated ground 
state the main part of the optical and EELS response in- 
tensity is determined by the excitonic sector, while for a 
band-like weakly correlated ground state it is the sector 
of unbound electron-hole pairs. Moreover, in the case of 
the strongly-correlated insulating cuprates it seems mis- 
leading to search for separate contributions of excitons 
and electron-hole continuum 1^ 

Concluding, it should be noted that the quantum- 
chemical CUO4 cluster model represents a physically clear 
albeit rather simplified approach to consider electron- 
hole excitations. However, it seems that such an ap- 



proach allows to catch the essential physics of the charge 
transfer transitions, and should be a necessary step both 
in the qualitative and semi-quantitative description of in- 
sulating cuprates. 
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